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Abstract 

The paper has three parts. In the first part we apply the theory of commuting pairs of (pseudo) difference 
operators to the (formal) asymptotics of orthogonal polynomials: using purely geometrical arguments 
we show heuristically that the asymptotics, for large degrees, of orthogonal polynomial with respect 
to varying weights is intimately related to certain spinor bundles on a hyperelliptic algebraic curve 
reproducing formulas appearing in the works of Deift ct al. on the subject. 

In the second part we show that given an arbitrary nodal hyperelliptic curve satisfying certain condi- 
tions of admissibility we can reconstruct a sequence of polynomials orthogonal with respect to semiclassical 
complex varying weights supported on several curves in the complex plane. The strong asymptotics of 
these polynomials will be shown to be given by the spinors introduced in the first part using a Ricmann- 
Hilbert analysis. 

In the third part we use Strebel theory of quadratic differentials and the procedure of welding to 
reconstruct arbitrary admissible hyperelliptic curves. As a result we can obtain orthogonal polynomials 
whose zeroes may become dense on a collection of Jordan arcs forming an arbitrary forest of trivalent 
loop- free trees. 
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1 Introduction and summary 



The present paper deals with the asymptotics of certain (pseudo-)orthogonal polynomials, its formal 
properties and connections with algebraic geometry. In order to explain the framework, let us recall the 
main results for ordinary orthogonal polynomials [14] in a simple exemplifying case. Let V{x) be an 



p n (x) be the (real) orthogonal polynomials (OP) for this measure. 

One of the main goals of modern asymptotic analysis is to describe their strong asymptotic as we let 
n — ► oo while N —> oo at the same rate. 

This problem has been brilliantly solved in [6, 7, 14, 15] using the associated Riemann-Hilbert problem 
(RHP); indeed the matrix 



solves a RHP with jumps on the real axis and rather simple asymptotics at z — oo [18]. More importantly 
the solution of this RHP characterizes the orthogonal polynomials; therefore if one could solve the RHP, 
then he/she would immediately have access to all information on the corresponding OP. 

Without entering now into the details it suffices to recall that the large parameter N and the degree n 
of the polynomials enter explicitly and in a simple way the Riemann-Hilbert data and the asymptotics at 
x — oo. Therefore all the "complication" of the asymptotics is controlled by the (highly transcendental) 
solution of the RHP. 

The technique developed in the last decade of the last millennium, called "nonlinear steepest descent" 
method, consists in transforming this RHP into a simpler, "asymptotic" one which differs from the exact 
one by small controllable errors, with suitable uniformly small bounds as N — > oo. 

The main character of the method is the so-called G-function (which we will discuss at length in 
the paper). Here we only indicate that it should satisfy some general properties which guarantee the 
amenability of the Dcift Zhou steepest descent method. 

The main logic of approach in most literature is then: 

• fix the potential V; 

• try to find an appropriate G-function for the given potential; 

• implement the steepest descent method 

In a certain sense this point of view mixes a "forward problem" (finding the asymptotics of orthogonal 
polynomials for a given potential) with an "inverse problem" (reconstructing the matrix solution of a 
monodromy/jump problem). 



even-degree real polynomial bounded from below and consider the Hilbert space L 2 (M,e NV ^dx). Let 



Y(x) : 



Pn(x) (f>n(x) 
p n -l{x) <j> n -i{x) 




(1-1) 
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Our approach is "purely inverse" ; namely we assign certain asymptotic data in the large iV-limit (which 
-of course- satisfy some consistency conditions) and remount the (class of) orthogonal polynomials and 
potentials whose asymptotics matches our given data. In a subsequent paper [1] it is shown that for any 
(polynomial) potential V (possibly complex) and any "Stokes' data" for the orthogonal polynomials, it 
is possible to construct an appropriate G-function. Together, the present paper and loc. cit. completely 
solve the problem of the large N asymptotics for these semiclassical orthogonal polynomials. 

1.1 Asymptotics of generalized pseudo-orthogonal polynomials 

In order to explain the main theorem of the asymptotic analysis we need to introduce some notations. 
The setting in which we move is rather algebro-geometric; the main piece of data is a hyperelliptic (nodal) 
curve C of genus g 

2g+2 

y 2 = M 2 (x) n (!"2) 

satisfying the conditions of Boutroux (we adopt the terminology of [24]) and admissibility (to be 
explained presently). Such curve is the asymptotic spectral curve; the reader acquainted with the 
asymptotic analysis of ordinary OPs may think of this curve as the hyperelliptic curve associated with 
the equilibrium measure. The Boutroux condition is a transcendental reality condition 



ydx e iR (1-3) 

'7 

where 7 is any closed loop on the spectral curve £. The Boutroux condition implies that (see Section 4) 
the set 

fio := (ft / ydx = 



is well-defined independently of the branch-point aj used in the integration; it also follows that it consists 
of a forest of (open) trivalent trees whose branches are Jordan arcs (possibly of infinite length) . 

It is then proved in Sec. 4 that this set uniquely defines a collection B of finite Jordan arcs joining 
the branchpoints ay which can be used as branch-cuts for the algebraic function y(x); one can then define 
uniquely (up to overall sign) the function 



h(x) :=U ( yx 

J ai 



in such a way that it is continuous on C and harmonic onC\S (see figure 1). 

The admissibility condition requires that the domains of negativity of (one of the two branches of) 
h be such that h < on both sides of each branch-cut (and also that the zeroes of M(x) do not belong 
to fi ). 
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An admissible Boutroux curve defines a potential via 

V'(x) = 2(y(x)) pol 

where the subscript means the polynomial part and the 
branch is such that h(x) ~ ^kV{x) at x = oo. 

1.1.1 Connection to the asymptotics of OPs 

Such an admissible Boutroux curve is connected to the 
asymptotics of generalized OPs for the potential V(x) above, 
in the following way (we omit some details for the sake of 
clarity; the full discussion is in Sec. 4). 

Given an admissible Boutroux curve, we can define some 
contours Tj satisfying the properties 

• They extend from oo to oo and approach it along di- 
rections for which ^RV(x) — > +00; 

• The (admissible branch of the) function h[x) is non- 
negative on these contours. 

In a figurative way each contour must extend from 00 to 00 
along different asymptotic directions and in order to do so 
is obliged to pass through the branchcuts B in order not to 
violate the second condition above. 

To each contour we associate an arbitrary 3 complex num- 
ber itj and then we can consider the generalized OPs as in 
[2] 





1 . 




/ / - 
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Figure 1: An example of admissible 
Boutroux curve, with the set fjo in evi- 
dence (solid-black contours). The branch- 
cuts defined by fjo are the thick arcs. The 
domains of negativity of h are the shaded 
regions. In this case the potential has de- 
gree 5 (see Sec. 5 for an explanation of 
why). The contours in texture represent 
choices of the contours Tj. In this case 
there are three inequivalent choices. For 
example the dashed contour comes from 
the left, follows two arms of the cross of 
thick arcs and goes to infinity downwards. 
The dotted lines are critical levelsets of h 
passing through a saddle-point. 



y^Xj / p n (x)p m (x)e NV{x) dx = 5 n 



(1-4) 



It should be clear that if we fix N, the contours Tj can actually be deformed quite arbitrarily in the 
finite region of the complex plane; it is only in the large iV-limit that we should require h < on these 
contours (this is the steepest descent condition). 

The main result of this part of the paper is Thm. 4.1 which provides strong asymptotic results for 
PN+r (where r is an arbitrary fixed integer) in terms of the solution of a model Riemann Hilbert problem 
(eqs. 4-21 in Lemma 4.6); This is the same logical approach used in [14, 15]. 

Such a RH problem is solved explicitly in terms of Theta functions and spinors on the Boutroux 
curve: this solution is essentially contained in [15] and is a special case of the more general approach in 



3 These numbers must be sufficiently generic. 
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[25]. However our approach is yet different and contained in Sec. 3 where we solve rather an inverse 
spectral problem (with a condition of Serre duality) and we find that the solution solves the desired 
RHP (see the next section in this introduction for the conceptual approach) . Note that the formulae that 
we obtain are significantly more streamlined compared to those appearing in the seminal work [15] and 
also "explain" the natural appearance of spinor bundles in terms of Serre-duality. 

The main point of this inverse asymptotic study is that Boutroux and admissibility conditions are 
sufficient (and also necessary as follows from the forward problem) for the amenability of the Dcift-Zhou 
steepest descent method 4 

A consequence and feature of this inverse asymptotic study is that we can construct very easily 
very complicated and somewhat "surprising" examples: we mention that a consequence of the uniform 
asymptotics given by the steepest-descent method is that the zeroes of the OPs condense on the Jordan 
arcs of the branch-cuts B in the complex plane; these arcs may form an arbitrarily complicated forest 
of loop-free trivalent trees 5 . 

The implementation of the steepest descent method requires a modification of the usual Airy parametrix 
near the turning points where three edges of a tree are incident (see Sec. 4.5). 

1.2 Boutroux and admissibility 

The Boutroux condition is a transcendental constraint on the coefficients of the polynomial defining the 
(nodal) hyperelliptic curve. At first sight it may not be clear that such Boutroux curves do exist and 
even less clear how constrained is the topology of f)o- 

This problem is addressed in Section 5; once more the approach is completely "backwards" . We start 
with an relatively brief summary of Strebel theory of quadratic differentials [32] adapted to our simple 
situation. This allows us to introduce Strebel's coordinates in which the Boutroux condition is just a 
linear constraint, hence trivially satisfiable. The importance of quadratic differentials for the study of 
2x2 Riemann-Hilbert problems has been pointed out in other papers (e.g. [27, 28]) but we are unaware 
of a systematic application to the steepest descent method as extensive as in the present paper. 

These Strebel coordinates parametrize different "cells" of different topology of Sjo but of the same 
dimension; in particular we can -given a topological forest graph satisfying some simple requirements- 
reconstruct by cutting-and-pasting a Riemann surface with the desired properties. This process is 
known in the literature of quadratic differentials as welding [32] : in our simple application the welding 
merges some half-planes and infinite strips into a Riemann-sphere. 

The condition of admissibility is also easily imposed in this approach since it simply singles out cells 
of specific topology. In a way in that section we are solving a more general problem of characterizing 

4 The nonlinear steepest descent method was suggested by Deift and Zhou in 1993 for oscillatory Riemann-Hilbert 
problems appearing in soliton theory. Some elements of the method appeared in earlier works of A. Its (1980-1982) [20] and 
also in [18]. An essential extension of the method, namely the introduction of the concept of G-function, allowed to include 
into the scheme in particular the RHP for orthogonal polynomials and it is due to Deift, Venakides and Zhou [14, 15]. 

5 Here "forest" is in the sense of graph-theory. 
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arbitrary polynomials by the topology of the graph of their critical horizontal trajectories (see Sec. 5 for 
the details). 

1.3 Spinor bundles and commuting difference operators 

In this first part we deal with formal and heuristic aspect of the asymptotic analysis. The formulas that we 
obtain will be shown a posteriori to represent the strong asymptotics of the orthogonal polynomials using 
the nonlinear steepest descent method. Of course the arguments used should be regarded as "temporary" , 
pending a rigorous proof (that is contained in the second part of the paper) . 

The main idea is the following: it is well known that any sequence of OPs satisfies a three-term 
recurrence relation which can be written in semi-infinite matrix form with a tridiagonal symmetric 
matrix which we denote by X. For generalized OPs w.r.t. complex moment functionals [2] the same is 
true but the matrix X has complex coefficients. 

If the potential is "scmiclassical" (i.e. if V'(x) is a rational function and the integrations supported 
on arbitrary arcs, see [2] for details) the quasi-polynomials ip n :— p n e~ N ^ satisfy a differential recurrence 
relation as well, with coefficients that form another matrix Y which has some closer relationship with X. 
In the simplest case of polynomial potentials V(x) then 

-^[^ ,^i,...] t = Y^o,^i,...] t (1-5) 
with Y antisymmetric and finite-band 

Y=\{V'{X) u -V'{X) t ) (1-6) 

the subscripts denoting the "upper" and "lower" part of the matrices. The commutation relations 
[jjd x ,x\ — jr imply the so-called string equation 

[X,Y] = il. (1-7) 

An heuristic approach to the asymptotics is that of sending N — > oo and also "look" at these matrices 
down the diagonal; if we shift the in dices of the matrices by N and send N to infinity then these matrices 
become doubly infinite. Moreover the string equation turns into a commutativity equation [Y,X] = 0, 
while retaining their band-structure and symmetry. 

In this regime these matrices can be regarded as commuting finite difference operators, a class 
of objects studied in several places in the literature. Indeed there exists a fairly general classification of 
commuting difference operators [26] which is the starting point of this part of the paper. 

In the case at hand the main data are: 

• a smooth hyperclliptic curve C of genus g with X being an invariant meromorphic function of degree 
2 with (symmetric) poles at two points; 
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• an arbitrary meromorphic function Y with prescribed pole structure 6 

• a suitably generic effective divisor T of degree g. 

These data allow us to reconstruct the matrices X and Y as the matrices expressing the multiplication 
operators by X and Y on a suitable basis of meromorphic functions. 

We extend this approach by introducing the Serre dual basis (of meromorphic forms) and a pairing 
(implementing Serra's duality) which takes the simple form of a residue. 

At this level the connection to OPs is still rather distant: in order to close the gap we need to 
consider a "self-dual" picture in which the meromorphic functions of [26] and our Serre-dual counterparts 
are replaced by sequences of spinors which are self-dual with respect to Serre-duality. This forces 
conditions on the divisor T and determines it up to the choice of a half integer characteristic. 

The actual application to the asymptotics of OP comes after tensoring this spinor bundle by a suitable 
flat line bundle; such bundle is provided automatically by the Ricmann Hilbcrt analysis and by the 
Boutroux condition. 

Indeed the Boutroux condition in this setting is § YdX G iM. for all closed loops; this means that 
exponentiation of J YdX provides naturally transition functions for a unitary line bundle with characters 
given by the periods of the Boutroux differential YdX. This last step provides the backbone of the 
asymptotic analysis of Sec. 4 

It should be added that this algebro-geometric part of the paper is more general and should yield the 
main ingredients for the asymptotic analysis of the most general semiclassical OPs, i.e. with potentials V 
whose derivative is rational and with integrals restricted to arbitrary arcs (hard— edges). However the 
technical details of the implementation of Deift-Zhou steepest descent method (Sec. 4) become rather 
more involved and we prefer to stay in the more standard setting of polynomial potentials for that part 
of the paper. 

Acknowledgements The authors would like to thank Peter Zograf for explaining the theory of Strebel 
differentials and Andrew Mclntyre for several discussion on topology of critical trajectories, Kenneth T.R. 
McLaughlin for several suggestions about relevant literature and Peter Miller for explaining how to obtain 
the density of zeroes of ordinary OP in the multicut case. Additionally we would like to thank the referee 
for suggesting several improvements to the exposition and historical background. 

2 Notation and main tools 

For a given smooth genus-g curve C with a fixed choice of symplectic homology basis of a and 6-cycles, 
we denote by ujj the normalized basis of holomorphic differentials 




(2-1) 



6 This is actually a generalization of the setting of [26]. 
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We will denote by 6 the theta function 



6(z) := ]T c 



i7rn-Tn—2iiTZ-n 



(2-2) 



The Abel map (with base-point p ) is defined by 



L 



u(p) = 



U) 



(2-3) 
(2-4) 



" rP rP 
/ Wl,..., / < 
-Jpo Jpa 

and is defined up to the period lattice Z + r • Z. For brevity we will omit any symbolic reference to the 
Abel map when it appears as argument of a Theta function: namely if p e C is a point and it appears as 
argument of a Theta-function, the Abel map (with a certain basepoint) will be understood. 
We denote by K the vector of Riemann constants (also depending on the choice of the basepoint) 



u) e (p) / u>j(q) - 5ji 



'33 



Pa 



(2-5) 



where in this expression the cycles a,j are realized as loops with basepoint po and the inner integration is 
done along a path lying in the canonical dissection of the surface along the chosen representatives of the 
basis in the homology of the curve. 

The Riemann constants have the crucial property that for a nonspecial divisor T of degree g, T = 
Sf=i7ji then the "function" 

/(p) = e(p-r-/c) (2-6) 

has zeroes precisely and only at p — jj, j = 1 . . . g. 

We will also have to use Theta functions with (complex) characteristics; for any two complex vectors 
e, 5 the theta function with these (half) characteristics is defined via 



e 



(z) := exp ( 2iir 



e-T-e 1 1 



e 



S e 
Z+ 2 +T 2 



Here the (half) characteristics of a point are defined by 

2z = 5 + re 



(2-7) 



(2-8) 



where the factor of 2 is purely conventional so that half integer characteristics have integer (half)- 
characteristics. In the sequel we will always use these half-characteristics. This modified Theta function 
has the following periodicity properties, for A,fieZ 9 



e 



(z + X + r/i) = exp iir(e ■ A — S ■ n) — iir[i ■ t ■ fi— 2inz ■ /i 







(z) 



(2-9) 
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Definition 2.1 The prime form E((, £') is the (—1/2, —1/2) bi- differential on L x C 

E(U')= ®^ C ~i\ (2-10) 



MC) 2 :=£a Ufc ine A 



w*(C) , (2-H) 

u=0 



fc=l 

where Uk are the normalized Abelian holomorphic differentials, u is the corresponding Abel map and 
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is a half-integer odd characteristic (the prime form does not depend on which one). 



The prime form E(£, £') is antisymmetric in the argument and it is a section of an appropriate line bundle, 
i.e. it is multiplicatively multivalued on C x C; indeed we have the multiplicative multivaluedness 

E(t+ aj ,C) = E(t,C) (2-12) 

£;(C + 6„C') = ^(C,C')exp^-^-^ C co^j (2-13) 

In our notation for the (half)-characteristics, the vectors a, (3 appearing in the definition of the prime 
form are actually integer valued. We also note for future reference that the half order differential h A is 
in fact also multivalued according to 

h A (p+a j )=e i ™*h A ( P ) (2-14) 
h A (p+bj) = e-^h A (p). . (2-15) 

Following the common usage in the literature, for a meromorphic function F (or section of some line- 
bundle) we will use the notation (F) for its divisor of zeroes/poles [16]. The writing 

(F)>-kp + mq (2-16) 

means that F has at most a pole of order k at p and a zero of multiplicity at least m at q. 

For C a hyperelliptic algebraic curve of genus g realized as a double cover of the plane, we will denote 
with p* the image of the point p € L under the hyperelliptic involution that interchanges the two sheets. 

To conclude this section we recall some expressions which are well known [17] and will be used later. 
Suppose that the hyperelliptic surface is given by the equation w 2 = Ilj=l 2 (' T ~ a j) an( l ^ us denote with 
oo± the two points above x = oo where w ~ ±x 9+1 . If we write the normalized first-kind differentials 
uij as 

Ex k ~ 1 dx ,„ _. 

°jfc ~ (2-17) 

k=l 



then it is a direct check to verify that 



e A (p-™+) = — (l + Oix- 1 )) , C A :=-J2<rj, 9 dj®A{0) (2-18) 
X j=i 
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3 Part I: Jacobi matrices and difference operators 



We follow [26] and consider the following setting: 

1. An hyperelliptic Riemann surface C of genus g 7 with X invariant under the holomorphic involution 
and with divisor (X) > — oo + — oo_. 

2. A meromorphic function Y with divisor 

(¥")>-/", fc^QOjn fc_ oo- + oo + (3-1) 

where the £j's are chosen in some subset of the Weierstrass points dX(^j) — O 7 . For later use we also 
point out that we could equivalently consider a meromorphic function such that Y(p*) = —Y(p); 
indeed we may "antisymmetrize" the above function by adding a rational function of X. This 
symmetry will be of use later on so we will introduce the special notation 

Y(p):=±(y(p*)-Y( P )) (3-2) 

3. A nonspecial divisor T of degree g, T = J2j=i 7r 

The genericity assumption on T is that all the divisors T + r(oo + — oo_) are also non-special 

i(r + r(oo+ -oo_)) = (3-3) 

where i(J)) denotes the dimension of the space of differentials with divisor exceeding D [16]. 
The genericity condition (3-3) implies also 

i(r+(r- l)oo+ - (r + l)oo_)) = 1 . (3-4) 

Indeed, if we had i(r+(r — l)oo+ — (r+l)oo_)) > 2 then, by Ricmann-Roch's theorem, there would exist a 
meromorphic function with (/) > — r — r(oo+ — oo_)+oo++oo_ . But then we would have that both / and 
X f (clearly linearly independent) would have divisor > — r— r(oo + — oo_) so that r(— r— r(oo + — oo_)) > 2 
and -again by Riemann-Roch's thm - we would have i(r + r(oo + — oo_)) = r(— r — ?*(oo + — oo_)) — 1 > 1, 
a contradiction with the assumption (3-3). 

By Ricmann-Roch's theorem and condition (3-3) it follows that for each r G Z there is a unique (up 
to multiplicative constant) meromorphic function with divisor 

(P r ) > -r - roo+ + roo_ (3-5) 

' Note that there is no loss of generality in assuming that X(oo M ) ^ X(ac u ) (i.e. the poles of Y have distinct X 
projection) because we may add to Y a suitable rational function of X to reduce ourselves to this situation. This is why 
we did not put any pole at oo+; moreover we can always add a constant to Y so that it has a zero at oo+. 
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The second condition (3-4) states that for each r € Z there is a unique (up to multiplicative constant) 
meromorphic differential F n with divisor satisfying 

(F r ) > r + (r- l)oo+ - (r+ l)oo_ (3-6) 

From these properties it follows that the product F r P s is a differential with at most two simple poles 
when r = s and otherwise it has only one pole. Therefore 

res P r F s oc S rs (3-7) 

oo_ 

This condition is just a manifestation of Serre duality. 

Since -for r — s- the product P r F s is a third kind differential with simple poles, its residues cannot 
vanish and hence we can normalize the two dual sequence so that their their Serre pairing (3-7) is actually 
o rs . 

Expressions for the functions P r and differentials F s in terms of Theta functions can be obtained 
following standard references, for example [26]. Since P r , F s play only a temporary role in our discussion, 
we delay explicit formulae until we arrive at the final objects of interest. 

The two sequences can be normalized in such a way that 

res F s P r = S rs . (3-8) 

oo_ 

Given the pole structure of X(p) we see that XP r is a linear combination of P r +i, P r , P r —i- indeed 

(XP r ) > -r- (r + l)oo+ + (r- l)oo_ (3-9) 

and the dimension (again generically) of the space of meromorphic functions with divisor exceeding the 
above one is r = 3 and it is spanned by the above three meromorphic functions [26]. 
We can express the coefficients of this three-term recurrence relation 

XP r = 7r+ iP r+ i + !3 r P r + 7 r P r _i (3-10) 

in two distinct ways: first and foremost, using the duality (3-7) 

jr+i — res F r+ iXP r , (3 r = res F r XP r , 7 r = res F r _\XP r . (3-11) 

oo_ oo_ oo_ 

This provides an explicit expression in terms of Theta functions if we express X as well in terms of them 

X( P ) = X n e f p - Z f n ^-^\ (3-12) 

where z ,zi are the two zeroes of X and X is a constant. A second independent way is obtained as 
follows 



XPjp) = - /r> Cr . r r det 

&et(P r+J (z k )) ]k=12 



Pr+lip) Pr(p) P r -l(p) 
Pr+l( z o) Pr( z o) Pr-l( z o) 
P r+ l(zi) P r (Zl) P r -i(zi) 



(3-13) 
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where zq,Z\ are the two zeroes of X (interchanged by the hypcrelliptic involution). The constant c„ is 
expressed in terms of Theta functions by matching the behaviors of both sides at one of the two infinities 
oo± and using (3-12). The dual sequence F n satisfies -by duality (3-7)- the transposed recurrence relation 

XF r = j r+1 F r+1 + r F r + 7r F r _ 1 . (3-14) 

3.1 Kernel and "Christoffel— Darboux" pairing 

We want to define the sequence of "projectors" K r r £ Z formally as the expressions 

r-1 
j=-oo 

so that res K r (p,£)P s (£) = P s (p) if s < r — 1 and zero otherwise. Clearly these expressions make little 

oo_ 

sense as they stand since they involve infinite series whose convergence should then be proved: what we 
want to have is a kernel K r (p, £) which is a differential in £ and a function in p that satisfies the following 
properties: 

1. As a differential in £ it has 

(a) a zero of order r — 1 at oo + and a pole of order r at oo_ 

(b) zeroes at T 

(c) simple pole at £ = p (the diagonal) with residue +1. 

2. As a function of p it has 

(a) A pole of order r — 1 at oo + and a zero of order r at oo_ 

(b) Poles at T 

(c) A simple pole at p = £. 

These properties define it uniquely as (see [17] pag 27. for similar kernels) and an explicit formula can 
also be written in terms of Theta- functions. Once more we will only write formulas for the final objects. 

The kernel K n (p,£) enjoys the Christoffel— Darboux property (an exercise using the recursion re- 
lations or the divisor properties) 

(X(p)-X(0)K r (p,0 = [F r -i(0,Fr(0][\ 7 o] [l^pf]^ ( 3 " 16 ) 

= 7r P r (p)F r _ 1 (0 - 7 r P r _i(p)F P (0 (3-17) 

Therefore we have the 
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Theorem 3.1 ("Christoffel Darboux" theorem) 

KrM = x{p) m (3-18) 

Taking the residue 

1 = res K r ( P , = 7rPr(p)J-r-l(p)-7rPr-l(p)fi(p) ^ 

i=p dX (p) 

we have a representation of the differential dX as 

dX = lr P r F r -i - 7 r P r _iF,. (3-20) 

3.2 Flat line bundles 

In the asymptotic analysis of orthogonal polynomials using the Riemann-Hilbert method that follows, 
we will need to tensor the line bundle whose sections correspond to the meromorphic functions P n and 
the Serre-dual line bundle (the F r 's) by a suitable flat line bundle. 

This line bundle can be described in much more general terms as associated to an arbitrary second- 
kind differential 8 . 

More precisely let rj be a meromorphic differential such that all residues are zero (second kind differ- 
ential) or integers. Let J r\ be its Abelian integral; it is defined on the universal covering of the curve C 
less the poles of 77. 

Definition 3.1 Near a pole c of n we define the exponential singular part E^^ of f rj as 

E V:C (p) = exp U n(0 ln(z(0 - z(p))j (3-21) 

where z(£) is a local coordinate z(c) = and the integral is along loop surrounding £ = c such that p is 
outside the loop. The function is independent of the choice of local parameter up to multiplication of a 
holomorphic function with nonzero value at c. 

Note that, since the residue of 77 at c is at most an integer, this "quantization" makes irrelevant which 
branch of the logarithm is used. 

The twisted line bundle £ r) . Associated to ry there is a line-bundle £ v with transition functions E n>c 
at the poles of 77 

Tensoring by £, v the line bundle described by T means that the sections of the line bundle will be 
"functions" ip with the following properties 

1. Poles at T 

2. Near a pole c of r\: ip n (p)E ri ^ c (p) — 0(1). 

The formulae for the dual sequences of wave functions/forms require minimal modifications and the 
expressions in terms of Theta function is an exercise that we delay to the next section. 

s We could in fact use arbitrary meromorphic differentials but this would introduce some slight additional complication 
in the formulae, and we leave this to another publication. 
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3.3 Spinors and the symmetric picture 

In the applications stemming from the (formal) asymptotics of orthogonal polynomials the sequence of 
OP should be orthogonal to itself: i.e. we should put in some way the two dual sequences on a symmetric 
footing; moreover the hyperclliptic involution should yield the dual sequence directly. This "symmetry" 
requirement fixes the divisor T. A similar construction in the theory of the algebro-geometric solutions 
of soliton equations goes back to the works [9, 10]. Indeed, let us look for a spinor s with the properties 

1. it has a simple pole at oo + 

2. it has simple zeroes at T 

3. no other poles or zeroes. 

Such a spinor exists provided that T satisfies some condition to be specified below (note that the degree 
of the above divisor is the correct one, g — 1): the square of a spinor must be a differential (with divisor 
of degree 2g — 2) that has divisor 

(s 2 ) = 2r - 2oo+ (3-22) 
Since the image in the Jacobian of the canonical class is —2K. we must have 

2r - 2oo + = -2JC , ^r + /C = oo + + ^ (3-23) 

2v = (3-24) 

[these equations are written understanding the Abel map and modulo the lattice of periods]. This 
determines the divisor T -up to the choice of a half-period v— from Jacobi's inversion theorem. We then 
define 

F r 

ir r :cx P r s , 7r* :cx — (3-25) 

5 

8 r (p,S)=8{p)K r (p,Z)^r (3-26) 

Exploiting the holomorphic equivalence between flat bundles and unitary flat bundles we can express 
s as section of the tensor product with the unitary line bundle characterized by 

X(a j ) = e i * A *, xfo) = e*"*' , (3-27) 

where A, B are the half-characteristics of v 

2u = A+tB. (3-28) 

Since v is a half period, then A, B £ I 9 . 
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3.3.1 O functional expressions 



In this paragraph we provide the explicit expressions of the spinors ir ni since they will be the final object 
of interest in our application. 

Both Tr r ,ir* are spinors (half integer differentials) belonging to a square-root of the canonical bundle 
of the curve tensored by the line-bundle £, v with transition functions E n _ c . and with divisors 



(n r ) > — (r + l)oo + + roo_ , (7r*) > roo + — (r + l)oo_ 

TT r (P+-f) = X(lUr(p) , 7T*0 + 7) =X~ 1 ("f)K(P) , V7G7Tl(£). 

with x(t) defined (in a basis) by eq. 3-27. 

Their explicit expressions are obtained using standard arguments and are given below 

i e^( P -oo_) 



/^e^ +1 ( P -oo_) 
i 

n 



A + e 
B + 6 

-A -I 
-B-S 
1 

" ~2 



(p + roo_ — (r + l)oo + )/iA(j>)e ^i' 
(p — (r + l)oo_ + roo + )/iA(p)e^ Q i ^ 



The constants h n are just the suitable normalizations so that 

res ir r ir* — 5 rs 

OO-j- 

They are 

6a(oo+-cx)_) 



e 



A+e 
B+S 



(r(c»_ — oo + ))0 



-A-e 
-B-S 



((r + l)(oo+-oo_)) 



and the spinorial kernel of the Christoffel-Darboux projector now reads 



@a(£ - oo + )6 A (p - oo_)~ 


r e 




(r(oo_ — oo + ) +p — £) 


6 A (C - oo_)6a(p- oo + )_ 




-4+e" 
B+(5_ 


(r(oo_ — oo + )) 



(3-29) 
(3-30) 



(3-31) 



(3-32) 



(3-33) 



(3-34) 



In the following applications 77 will be an antisymmetric differential (under the hyperelliptic involution); 
in this case we have the symmetry 

7T r ( P *) = -Tt*(p) (3-35) 

which explains the notation; indeed this follows from the fact that the Abel map (based at a Weier- 
strass point) of p* is the opposite of that of p and from the symmetries of the Theta functions with 
characteristics 



9 



A 
B 



(z) = e 



-A 
-B 



(-*) 



(3-36) 
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together with the antisymmetry of rj{p) = —T}{p*). With these normalizations the recurrence relations 
become automatically symmetric because 

7,- = res ir r Xir* +1 = — res ir*Xir r+ i — — res it*Xit r+ i = res Tt*X-n: r+ i — j r (3-37) 

00+ 00+* oo_ 00 + 

and hence 

Xn r = 'Jr'Kr+1 + /3 r TT r + 7 r _x7T r —i (3-38) 

Moreover it follows from the previous Christoffel-Darboux identities (Thm. 3.1) that 

7r r _i(p)7r*(f) - n r (p)K-i(€) 



dX{p) = J r ( 7T r _l(p)7T*(p) - 7T r (j?) 7T*_ 



(3-39) 
(3-40) 



3.4 Ladder matrices and Lax matrix 

We preliminary point out that any spinorial Baker-Akhiezer function irk can be written as a linear 
combination in terms of any other two consecutive wave-functions (forms) with polynomial coefficients 
in X: indeed the recurrence relations (3-38) can be rewritten in matrix form as 



n r+ i 



7I>+1 




7T r 





II r =: a(X)II r 



1 

The ladder matrix & r (X) is invertible and the inverse is linear as well in X. Therefore 

r-\-k 

Tl r+ k+i — a • • • all r =: a n r , k > 

r+k r r 

1 a - x n 



n r _ fc _i = a 

r-k-1 



r-1 



Even more directly, denoting by A r 



7r 
-7r 



k>0. 

the kernels K r are written as 
X(p) - X($ 



(3-41) 

(3-42) 
(3-43) 



Then we have the identity 

U {p) = res ^(O^feO = -E resH^O^^O = £ res "ff"*'^' 



n r (p) 



(3-44) 



(3-45) 



which follows from the fact that in any local parameter & r (p,0 = ^3i^32(l + 0{z — z'), with z = 
z(p), z' = z(£). It is easily seen that the matrix 



(3-46) 



is a polynomial in a; of degree |i — j\ and computes directly the product of the ladder matrices above. 
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3.4.1 Lax matrix 



Although Y II. r in general does not have the same pole structure as the sequence of the ir r 's, nevertheless 
we can express it in terms of the same vector II r . 

Let us define the rational function V' a of x, a = [M, — such that 



Y(p)-V^(X(p)) = 0(l) nearoo, 
Y(p) - VL(X(p)) = OiX- 1 ) near 



(3-47) 
(3-48) 



The function V!_ is a polynomial of degree fc_ whereas are polynomials in x _^ oc of degree and 
without constant coefficient. Define then the total potential 



V'{x) := J2 KV) ■ 



(3-49) 



This rational function of x has the property that Y(p) — V'(X(p)) is analytic near the points oo M and 
oo_. Then we find (using that Y has a zero at oo+) 9 



yn r (p)=resY(OlIr(0*r(P,O 



res yn r (0^-(p,0 

U,oo_ 



(3-50) 



V res ?(0n r (0^(p,0-^/es ?(0n r (0^(p,0 = 

E /es V'(X(0)IIr(0^r(P,0-Z)«8f(0nr(0<«r(P»0 = 

E res ^(x(0)n r (O n ;3 )A ?ff ~E/ e f y(Onr(0^-(p.O 
Q r^_ « =oo ° a (pj - a (^) ^ c=c 3 



V'(X( P )) 



* » £ = OO a if... 



E res f(0n r (0^(p,0 





V'(X(p)) 



n r (p) — > res 

V res fmn (0 g^jW(g) 
^res y(C)ii r (5) x(p )-*(0 





-*(0 


v'(x(0) 


- V"(X(p)) 



(3-51) 



*(p) " X{t) 



n r (e)n*'*(OAn r (p)- 



(3-52) 



Summarizing we have -in matrix form 

yn r = 



yn r = flW(x)n, 
o o 



V'(x) 

-E 



E 



res ^W^A, 
X-x 



(3-53) 



'Recall that represent those Weierstrass points (branchpoints for X) which coincide with a simple pole of Y (or K). 
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Quite clearly, the matrix Dj 1 (x) is a rational matrix with the same pole structure as the derivative of 
the "potential" V'(x): at the other poles (£ r ) of Y, it has simple poles with a nilpotent residue. 

Similarly we could repeat the computation for the antisymmetric function Y; however it is immediate 
to see that 10 

Y = hf'{X) -Y (3-54) 

so that the matrix Di v \x) representing the multiplication by Y differs from the matrix Dr(X) only 
by the multiple of the identity \V(X)\ 

D to( x ) = YMi _ d^(x) , [U r ,U* r ]Ya 3 = D^(X)[U r ,U*} . (3-55) 

where we have used that II* solves the eigenvector equation with eigenvalue — Y since it corresponds to 
evaluation on the other sheet of the X-projection and Y is antisymmetric. 



The WKB Ansatz. We now spend some words on the heuristics of the whole construction and on 
its logical and historical underpinnings. As pointed out in the introduction (at least for polynomial 
potentials) the OPs' satisfy a three term-recurrence relation of the same form as (3-38) and also a 
differential recurrence relation (1-5). The compatibility relation between the two is contained in eq. (1- 
7), often referred to as the "string equation" for the recurrence coefficients [31]. It can be rewritten in 
matrix form as the compatibility relation between the following differential/difference equations 



N 



d x ^ n = D n (x)^ n ,* n+1 = 



x-/3„ 
1 



7n-l 




Pn ( X ) <f>n 



rV(x)<7 3 



(3-56) 



where n = N + r (with r 6 Z bounded as N — > oo), <j> n are the Cauchy transforms of e~ NV p n 's and 7 n ,/3„ 
denote the recurrence coefficients for the orthogonal polynomials (see for example [3] for explicit formulas 
for the matrix D n in a very general setting of semiclassical orthogonal polynomials). The heuristic 
approach is then to assume the validity of Wcntzel- Kramers- Brillouin (WKB) ansatz for ^„ = ^N+r m 
the form 

V N+r ~~ r e~ N ^r Y ^ (3-57) 

where S r (see eq. (3-68)) should be expanded (as well as Y) in formal series of A -1 . Plugging the Ansatz 
(3-57) into the differential equation (3-56) yields, to leading order 



(3-58) 



which is nothing more than an eigenvalue/eigenvector equation precisely of the type appearing in (3-55) 
(with 77 = NYdX). In short the WKB approach amounts to substituting the ODE in (3-56) by its 
symbol, thus turning it into a spectral problem. 



10 Indeed Y(p*) + Y(p) is invariant under the hyperelliptic involution and hence it is a rational function of X, namely 
V'(X) in our notation. 
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This line of approach was developed in the works [30, 31], based in turn on the isomonodromic 
techniques in the theory of Painleve equations [21]. In the context of OPs it was applied to the "one-cut" 
case in [6] and to multi-cut settings in [8]. 

The self-consistency of the WKB method calls for the periods of the exponent J x Y(x)dx to have 
no real part, for otherwise the asymptotic expansion in inverse powers of N for S_/v would not be valid 
under analytical continuation; such condition (Boutroux condition) will be the main requirement in 
the second part of the paper. 

In particular, as anticipated in the introduction, the asymptotic spectral curve (1-2) with the Boutroux 
condition (1-3) should be thought of as the "limiting spectral curve" of the matrices D n (x) appearing in 
the differential equation of eqs. (3-56). 

The matrix S r consists of genuine Baker-Akhiezer functions. The relation with spinors and the 
Serre duality (3-32) should be understood as a large N relic of the orthogonality relations: indeed, at 
a very formal/suggestive level, the orthogonality reads J p n p m e~ NV ( x >dx = S nm which suggests that a 
symmetric way of considering the large- N limit should be by formally splitting dx = VdxVdx so that 
PN+ r e~^ v ( x ^Vdx ~ 7r r . The iV-dependence will be in the asymptotic line bundle determined by the 
differential rj = NYdx on the hypcrclliptic spectral curve, while the finite perturbation r will give the 
lattice behavior (of Toda type) . 

Of course, the reader should regard these considerations as a mere back-of-the-mind motivation to- 
wards an educated guess, the full justification of which comes only a posteriori via a Dcift Zhou nonlinear 
steepest-descent analysis (contained in the second part of the paper) adapted to particular hyperelliptic 
curves satisfying the Boutroux condition. The existence/construction of these is the content of the third 
part of the paper. 

Once the strict connection between the 7r r 's and the uniform asymptotics over compact sets of the 
orthogonal polynomials are made rigorous, we could conclude that indeed det(yl — D n (x)) admits a 
bona-fide limit and this limit is a Boutroux curve. 

Indeed it will be shown in Section 4 that the solution ^N+r of the Ricmann- Hilbcrt problem for the 
orthogonal polynomials (eq. (3-56) is indeed approximated by the matrix z. r (x)e N f Ydxa3 (with S r given 
in (3-68)) uniformly over compact sets not intersecting the branchcuts. 

Thus we have that the matrix of the differential equation D N+r (x) = ^^'N+ r ( x ) 1 ^N 1 +r( x ) (which is a 
polynomial in x of constant degree equal to deg V^'(a;) for polynomial potentials) also tends uniformly to 
D r }{x) — 'R{x)~ 1 Y{x)a^B{x), which is our eigenvalue-eigenvector equation (3-55) for the line-bundle n — 
NYdX. Therefore a fortiori also the spectral curve must tend to the asymptotic Boutroux hyperelliptic 
curve. Since the convergence is uniform over compact sets and the objects D^+r and its characteristic 
polynomials are all polynomials in x of bounded degree, the convergence is established in any bounded 
set of the plane. 
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3.5 The spinor Ha 

Following ([17], page. 13) we can explicitly write the half order differential (up to a multiplicative 
constant which is irrelevent in our application since it would then be reabsorbed in the normalizations). 
Let us denote by ctj = X(£j), j = 1, . . . 2g + 2 the critical values of X (at the Weierstrass points) and 

2g+2 

W 2 := Yl ( x ~ a j) O 59 ) 

3=1 

There are 4 9 half periods in the Jacobian: they are all in one-to-one correspondence with the 4 9 partitions 
{ii, . . . ,i g+ i-2m} U {ji, ■ ■ ■ , j g +i+2,n} of the set {1,2, . . . , 2g + 2} for m > as follows: 

1. All even non-singular half-periods are given by the choice m = and are the image in the Jacobian 
of 

9+1 

cti k — oo + — oo_ — K, (3-60) 

fc=i 

2. all odd non-singular half-periods are given by the choice m — 1 and are the image in the Jacobian 
of 

9-1 

a y - K ( 3 " 61 ) 

3=1 

3. All other half-periods are singular and they are even or odd according to the parity of m. 

Any choice of odd, nonsingular half-integer characteristics corresponds to a choice of g — 1 points 
amongst the Weierstrass points, £i t ,... , £j 15 therefore 

h\^\[{X-a lk ) w (3-62) 

fc=i 

In other words, with proper understanding of the analytic continuation on the (double cover of the) 
hyperelliptic curve, the spinor h& can be expressed as (up to overall constant) 

h A = ^ 9k =^ X ~ aik) l VdX (3-63) 

Uk=i( x - a jJ 4 

3.6 Riemann Hilbert problem 

The spinors (3-31) exhibit a multiplicative behavior which depends on the half period v = + |r • B 

TT r (p + a j ) = e"- 4 ' ir r (p) , ir r (p + b ) = e l7lB > ir r (p) (3-64) 
n*(p + aj) = e-" A ^*{p) , <(p+ 6j) - e-^Tt* (p) (3-65) 

(3-66) 
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These phases define a character 



and in general so does any half-period (or half-integer characteristics). 
The matrix 



'dX 



7T r _l lTT*_l 



satisfies 



for Di v) (X) given by the formula (3-55). Moreover, from (3-40) it follows 



YE r a 3 = £>("> (X)Er 



det ^ r 



i 

It 



(3-67) 

(3-68) 
(3-69) 

(3-70) 



The matrix 3 r is a bona-fide multivalued matrix- valued function on the (desingularization of the) spectral 
curve 

= H(X, Y) = det(Fl - D^\X)) . (3-71) 

It can be thought of as a multivalued function of x = X(p) with branchpoint singularities at the branch- 
points otj = X(£j). 

On the sheet of the map x — X containing the point oo + it solves a certain Riemann-Hilbert problem 
which we describe below; the main complication arises from the proper understanding of the square-root 
of the differential dX. We note that the entries of H„ are proportional to (on each sheet) 

Clearly this "function" makes sense only after a suitable surgery on the plane, which is what we describe 
in the paragraph below; moreover it depends on the choice of the characteristic A under the identification 
between half-periods and partitions of the Weierstrass point described earlier. 
Each entry of S is of the form (refer to eq. 3-31) 



Qa(x) 



JdX 



(3-72) 



TT r 1 Q A (p — 00_) 

7dT = 7^Q a +1 (p-°o+) 

= : F r (p)Q A (X(p)) 

K , = i Qa(p~°°+) 

fdX" VKG r A +1 (p-oo- 
= : F*(p)Q A (X(p)) 



e 



e 



A+e 
B+6 



-A-e 
-B-8 



(p + roo_ -(r+l)oo+)e J «i n Q A (X{p)) 
(p-(r+ l)oo_ + roo + )e^i v Q A (X(p)) 



(3-73) 



where the functions F r and F* have the monodromy 

F r (p + 1 )=Xu+A(l)F r (p) 

f;(p+i) = x»+a(i)f;( p ) 



V7e7n(£) 



(3-74) 
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where we have used x„xa = Xv+A and X— A = XA (this last valid for half-periods only). It follows 
also from the above formulas that each entry of S r has a singularity at the branchpoints otj of type 
[x — fx,) -1 / 4 . In fact Oa(p — oo_) (and Oa(p — °°+)) nas simple zeroes at the Weierstrass points 
cti k = X(£i k ) ,k= 1, ... g — 1 appearing in eq. (3-72); since the local coordinate is y/x — oti k we see in 
eqs. 3-73 that at all branchpoints we have the advocated behavior. 

When thinking of F r ,F* as functions on a simply connected domain of C (one sheet of the X- 
projection), they define functions with singularities only at the branch-points and essential singularities 
at the X-projections of the poles of the twisting differential r\. 

If r\ is an antisymmetric differential (w.r.t. the hyperelliptic involution) 11 then there is the further 
symmetry (which follows directly from the explicit formula defining them) 

F r (x) = -F*(x*) . (3-75) 

(Nonstandard) surgery. Let B be a set of branchcuts for the projection X: specifically, if on are the 
critical values of X then the smooth hyperelliptic curve is written as 

29+2 

W2 = II ( X ~ ^) • ( 3 " 76 ) 

i=l 

Then B is a collection of mutually non- intersecting oriented arcs £j, joining two points aj 1 ,aj 2 in such 
a way that one branch of W can be defined as a single valued function on C \ B. The "standard" way of 
performing these cuts is to join otn to Q2i+i (whatever numbering has been chosen). We point out that 
there are "nonstandard" ways of performing an equally satisfying surgery; the only condition is that at 
each point cti originates an odd number of cuts (we will need this generality in the following). We also 
require (which will be enough for our later application) that C \ B is connected and that each connected 
component of B has an even number of vertices. We add some oriented arcs (called gaps) joining each 
connected component of B to the next, and the last one to oo. We denote by £ the collection of all 
oriented cuts B and oriented gaps (see for example Fig. 4). On the resulting simply connected domain 
C\Ewe have the following Ricmann-Hilbert problem. 
Near the poles c we have 

Sr(p) ~ near c (3-77) 

cr 3 := diag(-l, 1), X(p) = x. 

The jumps discontinuities of this Riemann-Hilbert problem are given in the following paragraph. 

11 There is not much loss in generality in assuming that it is antisymmetric, because it can be always antisymmetrized 
by an exact differential which does not change its characteristics. 
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Jumps on S. For x € £ \ {ai}i=i...2 g let 7(2;) be a closed positively oriented loop in C \ E intersecting 
sigma only at x (the simple connectivity of C \ £ implies that the homotopy class is unique) . 

Define $/\(x) as the difference between the number of Weierstrass points (in the interior region cut by 
7(2;)) entering in the numerator of Qa (which define A) and the number of the ones in the denominator. 
It has the following properties which are easily proved: 

• it is odd on the cuts; 

• it is even on the gaps; 

• if x is on the last gap extending to infinity, then Ua(x) = —4 = (g — 1) — (g + 3). 

The jumps of Qa(x) on the cuts and the gaps are given by e 2 ' r *( 1± iltA(a:)) depending on the orientations 
of the corresponding cut/gap. In particular the function Qa(x) is continuous across the last gap. 

Any Z2 (which we think of as the multiplicative group consists of 1 and —1) character on tti(jC) 
induces an assignment of signs on the cuts and gaps of £ as follows. Define a loop 7(2;) such that: 

• 7(2;) is a closed loop on C based at £1 (A(£i) = ot\) the Weierstrass point chosen as base-point for 
the Abel map; 

• the A-projection of 7(2;) is a positively oriented loop intersecting £ only at Q\ = A(£i) and 
at x. 

The above recipe defines 7(2;) on C up to orientation; however this ambiguity is inessential for us 
because we will be evaluating Z 2 characters only (and (— 1) _1 = —1). 

Given an arbitrary H : 7Ti (£) — > Z 2 we assign to each oriented cut Sj and each oriented gap T,g the 
sign 

(-fj := H( 7 (x))(A( 7 (2;)) • £;) , for x e S j 

(-fi := H( 7 (x)) for x E S £ (3-78) 

where (A (7(2;)) • Sj) denotes the intersection number for oriented curves in C. Some simple properties 
are worth pointing out: 

• If we denote by £1 the cut attached to ai, then for x <E Si the loops 7(2;) are homotopically trivial 
(in £), hence the character ^(7(2;)) = 1; 

• there are 2g + 1 amongst cuts and finite gaps (i.e. excluding the gap that extends to infinity); 

• as x moves along the links of E then 7(2;) spans 2g homologically independent loops (defined up to 
orientation) in C; 
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• any Z2 character on ni(C) corresponds (in one-to-one fashion) to a half-integer characteristic on 
the Jacobian of C. 



These facts imply that we can arbitrarily assign signs to all finite links of S (except Si) and then define 
a Z2 character on 7Ti(£) by using eqs. (3-78). This is consistent since the rank of it\(C) is 2g and there 
are (excluding Si) the same number of finite links in S. 

With this preparatory material we can formulate the jump relations of 5 on E 

Lemma 3.1 The matrix S r (a;) on C \ S satisfies the following jump relations: on the cuts 

E r (x) + = (-^^^XA+AXhix)) ■ ^)E r (x). ( _° : I ) (3-79) 

while on the gaps 

E r (x) + = (-l)^ Xl , +A ( 7 (x))~ r (:r)_ (3-80) 

The proof is a simple inspection of the properties of F r (x), F*(x) and Qa (using the argument principle). 
The main point that we make here is that on the gaps the jump is at most a sign and on the cuts it is 
given by ( ^ J j up to a sign (which is computed in detail by the formula). 

We have remarked above that we can assign arbitrary signs to all finite links of S except Si and lift 
this assignment to a Z 2 character of 7Ti (£) . Therefore we have proved 

Proposition 3.1 If we choose the orientation of Si such that (— 1) A 2 (X(j(x)) ■ Si) = 1 then, for 
any choice of A and of orientations of the other cuts/gaps, there is a corresponding (unique) choice of 
the half-period v such that the matrix H„ solves the sign— normalized Riemann-Hilbert problem 

on the cuts (3-81) 

(3-82) 




3 ' 



(3-83) 



namely it has no jumps on the gaps. 

4 Part II: asymptotics of orthogonal polynomials 

The construction in section 3.3 suggests that the spinors could be interpreted as the large N asymptotics 
of some orthogonal polynomials. In this section we will show that this is indeed the case. We will make 
use of the steepest decent method [11], [14], [15], [33], [19] to approximate the Riemann-Hilbert problem 
satisfied by orthogonal polynomials with semi-classical potentials [2] for large N. The solutions of these 
Riemann-Hilbert problems then represent the strong asymptotics of the orthogonal polynomials in the 
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large N limit. The main result of this section is that, when the meromorphic functions Y and C and the 
Riemann surface satisfy a certain condition (Definition 4.2), the spinorial matrix constructed in Sec. 3.2 
provides the solution to this deformed Riemann- Hilbcrt problem and hence gives the strong asymptotics 
of the corresponding orthogonal polynomials away from the branch points of Y where Airy asymptotics 
must be used instead. 

The setting is not as general as in the first part of the paper: the potential will be just polynomial, 
namely the divisor of poles of y = Y(p) will coincide with that of x — X(p). Additional complications 
would arise in presence of other singularities and hard-edges but very interesting and somewhat surprising 
features are already present for this simplest class of potentials: the general situation is addressed in a 
future publication. 

Let (£,x,y) be a triple such that £ is a smooth hyperelliptic surface of genus g, x = X(p) is an 
involution-invariant meromorphic function with two simple poles and y is a meromorphic functions on 
the Riemann surface C represented as a (nodal) hyperelliptic relation between x and y 

y 2 = U^+ 2 (x~ ai )M 2 (x) 

where the a* are the (distinct) x-values of the Weierstrass points and M(x) is a polynomial with roots 
not coinciding with any of the ay's. 

Let V(x) = ^f x% be a polynomial of order d+ 1 (without constant term) such that 

Y(p)~±V'(X(p)) + 0(l), p^oo+ (4-1) 
on C. The degrees of V, M and the genus are related by 

g + 1 + deg(M) = deg(V) . (4-2) 
We define the Stokes' ray for this triple (C,x,y) as follows 

(2k - 1)tt' 



7L k = reC,arg(j;)eiJ + 



2(d + l) 



*: = arS i" rf+l) , fc=0,...,2d + l (4-3) 

The situation which will be relevant for our discussion is the following, which we formalize in a 
definition 

Definition 4.1 The triple (C,x,y) is said to satisfy the Boutroux condition if all contour integrals 
§ ydx € tR are purely imaginary. 
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Since y is anti-symmetric under the hypcrclliptic involution and has only poles above x — oo, there 
are 2g + 1 homologically independent classes to consider. It will be proved later in Section 5 that such 
condition can be fulfilled. In fact these "Boutroux" triples will be constructed. 

Let ai be one of the branchpoints: we can define the following function on the curve C 



hip) = n ( r ; 

\J ax 



y(s)dsj . (4-4) 

The choice of a± or any other branchpoint does not affect the definition of h (see second bulleted item after 
Def. 4.2). Here we think of h(p) as a function on the hyperelliptic curve itself: the contour of integration 
is immaterial because of the Boutroux condition (the only additive monodromy of the integral is purely 
imaginary). Therefore h(j>) is a well-defined harmonic function on C \ {oo-t}. If we consider it as a 
function on the x-plane, then it has two branches which differ only by a sign: in particular the harmonic 
continuation of one branch on the punctured x-plane C\{aj}_ 7 - = i ) ... ) 2 9 +2 has the exact same multiplicative 
monodromy of the analytic continuation of y (again because of the Boutroux condition) . 

Lemma 4.1 Under the B outroux- condition, the function h(x) = !ft ( y(s)ds^ has only multiplicative 
monodromy with values ±1 and is otherwise independent of the choice of contour of integration where 
yi(x) is the branch of y that behaves like \V'{x) near x — oo. Its zero level set is well defined. 

The multivaluedness of h{x) is the same as the multivaluedness of y{x): therefore we can make appropriate 
cuts on the x-plane for which y(x) becomes single valued on the resulting domain. On the same domain 
then h will be harmonic, possibly with jump-discontinuities across those cuts. If x belongs to one of these 
cuts and h{x)± denote the boundary values on the two sides, we have 

h(x)+ = -%)_ . (4-5) 

If we can choose the cuts within the zero-level set of h, then h will be continuous on the whole plane 
and harmonic away from the cuts. For the time being we formalize this into the following definition. 

Definition 4.2 Let (C,x,y) be a triple such that C is represented as a nodal hyperelliptic curve by 

y 2 =M 2 (x)U 2 l+ 2 (x- ai ) (4-6) 

where ai are distinct and that M(x) is a polynomial with roots not coinciding with any of the otj 's. This 
triple is called admissible ( and noncritical ) if 

1. § ydx is imaginary for any closed curve 7 S C and res ydx — 1 (normalized Boutroux 

7 ' p— 00+ 

condition 12 ^. 



12 The adjective "normalized" refers to the normalization of the residue at 00+; clearly, given a non-normalized Boutroux 
triple one can get a normalized one by rescaling y by a real constant. 
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2. It is possible to define branch cuts of y in such a way that 

(a) all the cuts are finite Jordan arcs denoted by Ej joining two branchpoints; 

(b) h(x) is continuous on the whole plane and harmonic away from the cuts; 

(c) for all cuts h(x) (or —h) is negative on both sides of each cut. 

3. For fa any root of M(x) then ^ (honcriticalityj. 
Since the zero levelset 

£ := {x\ h{x) = 0} 

is independent of the harmonic continuation it makes sense to study its topological properties. 

Before proceeding, let us make a few observations of the set S)o in the presence of the Boutroux 
condition. 

• The set f)o consists of a finite union of Jordan arcs. 

• All branch-points ctj belong to fio'- indeed J" J ydx is half of a closed loop on the Riemann-surface 
C and -by the Boutroux condition- it is thus purely imaginary. Therefore h(atj) — 3? ydx = 0. 

• Since we assume all the branch points to be simple, there are 3 arcs originating from each branch 
point (as a simple computation in a local coordinate shows). They can either lead to another branch 
point or towards oo. We call a branch that ends at another branch point a closed arc (or branch) 
and one that leads to oo an open arc. Moreover, the only self-intersection points of the set f)o are 
the branch points. 

• The set Sjo cannot contain a closed finite loop 7: if this happened, then there would necessarily be 
another closed loop within the region bounded by 7 and without any branch point or singularity 
inside it. Since then (one branch of) h would be harmonic on this simply connected domain, 
continuous on its closure and with zero boundary value, h(x) would vanish identically by the 
maximum modulus theorem, a contradiction. 

• The function h(x) is continuous across a cut Ej if and only if Ej C fio- 

• There can only be an odd number of branch cuts coming out of a branch point for otherwise the 
continuation of y(x) in a neighborhood of the branch point would have no multivalucdncss and the 
point could not be a branch-point. In our situation we can run either one or three cuts at each 
branch-points while remaining in £) . If there is only one branch-cut then we cannot choose the 
branch of h so as to have a definite sign in a neighborhood of the branch-point. If there are 3 
branch-cuts, then we can choose a branch of h such that h is continuous (not harmonic) with a 
semi-definite sign in a neighborhood of the branch-point. 
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These simple observations imply that f)o is a trivalent graph with no closed loop with trivalent vertices 
at the branchpoints. It may contain possibly some open Jordan arcs not containing any branch-point 
(i.e. extending from oo to oo). We will dwell at length on the topology of such graphs in Sec. 5. 

The second condition in the definition of admissibility (Def 4.2) implies that all branch cuts belong 
to fjo (because otherwise h would not be continuous across the cut): therefore none of the branch points 
between different connected components of S)o can be connected by a branch cut. We then have 

Lemma 4.2 If(C,x,y) is admissible, then each connected component of Sjo contains an even number of 
branch points. 

Proof. We can form another graph B using the branch points as vertices and the branch cuts as 
edges. Let the connected components of this graph be Bi. We will show that each Bi contains an even 
number of vertices. 

First note that by the last bulleted item in the above list of facts, each branch point a, € Bi can only 
be connected to an odd number (specifically 1 or 3) of branch points through branch cuts coming out of 
it. Let [3 be a branch point connected to on through a branch cut: then f3 can only be connected to an 
odd number of branch points. Apart from on, all of these points can only be connected to cti through 
this branch point /?, or a closed loop would form in the graph. Therefore the branch points that are 
connected to (3 will add an even number of new branch points to the connected component of Bi of the 
graph. 

By repeating this argument, we see that the total number of branch points in this component must 
be even. Q.E.D. 

Definition 4.3 We will call a branch point with 0, 1,2 open arcs incident to it, a branch point (vertex) 
of type I, II or III respectively. 

In the set fjo there are open arcs which do not contain any branch-points: removing those arcs yields a 
trivalent graph Xo (the critical graph) with some open edges (i.e. edges attached to only one vertex). 
We can now define the branch cut structure of such a trivalent tree-like (or forest-like) graph. 

Definition 4.4 Let Xo be a connected trivalent graph with no closed loop and open end-edges. Then the 
branch cut structure of the graph Xo is a subgraph B of Xq containing all its vertices and such that 

1. Each vertex has either 1 or 3 edges incident to it, 

2. B has no open-edges. 

3. Each connected component of B has an even number of vertices. 

The edges of B will be called branch cuts. The definition applies also to graphs Xo with several connected 
components provided each component has the aforementioned properties. 
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Figure 2: A trivalent open graph with the properties of Def. 4.4 and its branchcut structure (thick lines). 

The following lemma shows that a branch cut structure for a graph with an even number of vertices 
exists and is unique. 

Lemma 4.3 Let Xq be a graph with an even number of vertices a%, . . . , ct2k as in definition 4-4- Then 
the branch cut structure B of Xq exists and is unique. 

Proof. We can reason independently on each connected component of Xq, so there is no loss of generality 
in assuming Xo to be connected. 

First remove all open edges and call Xo the result of this first pruning: it is still connected and with 
only trivalent, univalent and bivalent vertices. 

If the resulting graph Xo does not contain any bivalent vertices then B = Xq. 

Viceversa, if Xo contains some bivalent vertices, then each of these belongs to a maximal chain of 
bivalent vertices; the chain connects two subgraphs of Xo which we denote by L and R (also with only 
uni/bi/trivalent vertices). They connect to the chain by an oddvalent vertex (if the vertex of -say- L is 
univalent, this means that L itself is a single vertex). 

Now, if the chain has an even number of vertices, then either both L and R have an even number of 
vertices or they both have an odd number of vertices; we need to remove every second link in the chain 
so that each bivalent vertex is turned into an univalent one, and there is only one way of doing so such 
that all the resulting graphs contain an even number of vertices. 

If the chain has an odd number of vertices, then either L or (exclusively) R have also an odd number 
of vertices; once more there is only one way of removing every second link in the chain so as to leave each 
subgraph with an even number of vertices. 

The result of this second pruning is a collection of elementary segments (two univalent vertices joined 
by a link, which is what remains of the bivalent chain) and the two subgraphs L and R (possibly with 
one extra vertex on one of the two if the chain was odd). By repeating this procedure we could reduce 
the graph into a collection of subgraphs {X^} which do not contain any bivalent vertex. The branch cut 
structure B of Xq is then given by B = UjXj. Q.E.D. 
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In order to apply this construction to an admissible triple (C,x,y) we adopt the following strategy: 

• Define h(x) — f yids, where y\ is the branch of y that behaves like ^V'(x) near x = oo with 
some arbitrarily chosen branch cuts. We can then define ,f)o = {h(x) = 0}; 

• remove the arcs of ,f)o that contain no branch-point of y and call Xq the resulting trivalent tree-like 
graph. 

• The branch-cut structure B of % can be used to redefine the branch-cuts of y\ (whence the name 
of "branch-cut structure"). These branch-cuts may be "nonstandard". Indeed it may happen that 
a branchpoint has three branch-cuts connecting it to as many other branch-points. Although 
unusual this is consistent with the multi-valuedness of y. 

The admissibility requirements then is the condition that the locally harmonic function h could be chosen 
such that 

1. h is continuous on the x plane and harmonic away from the branch-cuts; 

2. the sign on both sides of each branch-cut is negative. 

At this point the reader may wonder how strong a condition is this and if there are any examples; in 
Sec. 5 we will show how to reconstruct an arbitrary admissible triple from the topological structure of 
its (admissible) graph Xo- For the time being we present only two figures exemplifying two curves of the 
same genus satisfying the Boutroux condition but one being admissible and the other non-admissible. 

Example 4.1 The first example is an admissible triple (left in Fig. 3, with 

y 2 = {x + l + r){x- 1 - r)(x - e i7r/3 )(x - e 2i7r/3 ) , ,r~ 0.4144... (4-7) 

The second curve is not admissible, although it satisfies the Boutroux condition (rather approximately, on the left 
of Fig. 3) 

y 2 ~ (x ~ 2 - 0.15i)(x + 2 - 0M)(x - l)(x - 0.88i) (4-8) 

For an admissible triple, we will choose branch cuts of y that satisfy condition 2 in Def. 4.2, which 
are the ones defined by the branch cut structure B of the graph Xo- 

4.1 The G- function 

The main character of the construction is (as in [19, 14]) the so-called G-function. In this context it is 
simply the Abelian integral J ydx; of course attention must be paid to its multivaluedness (which is only 
additive on the Riemann surface C). 

The construction of the previous section has yielded a forest graph B = UiBi made of several connected 
components each of which is a loop- free treelike closed graph with only odd-valent (1 or 3) vertices. 
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Figure 3: The critical graph Sj = X of two Boutroux curves: the one on the left is admissible (marked 
are the signs of h), the right is not admissible. 




Let us order these components in some arbitrary fixed way; using some arcs not intersecting B we 
then connect a vertex of B\ with a vertex of B2, then one of B2 to one of B3 and so on and so forth, 
adding an open arc that extends to 00 from the last connected component of B : we call E the resulting 
set (see Fig. 4). We will call these arcs joining the consecutive Bi's the complementary arcs or gaps 
and denote them by Ey, whereas we denote by Ej the edges of B. 

The set C := C \ E is simply connected and hence we can define unambiguously 

G«\x)= [ X yi(s)ds, i= 1,2, (4-9) 

where y\ — —yi are the two branches of y and ct\ is a univalent vertex of B\ and y\ ~ |V' (x) at x = 00. 
The contour of integration is taken to lie within the simply connected domain C. 

We also assume that the edges of B have been oriented, so as to be able to distinguish the left from 
the right side. We then formulate 

Lemma 4.4 The function G(x) :— G (x, E) has jump discontinuities on E as follows. 

G(x) + = G(x)-+i<Tj, x G Ej 

G(x)+ = -G(x)- + i(Tj, xgEj (4-10) 

where ai and ai are real. 

Proof. Let us first consider a gap Ej and x G Ej: the two paths from ct\ to x on the left/right of the 
gap lift to a closed loop 7 on the curve C because it encircles some number of connected components 
B\, . . . ,Bj, each of which contains an even number of branchpoints. Therefore 



G(x)+ = G{x)- + / ydx , (4-11) 
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which proves the first identity, with aj = —i J ydx. 

In order to prove the second identity we first note that 

G 1 (x) ± = -G 2 (x) ± on a branchcut. (4-12) 

Let x G Ej (an edge of B); a closed loop on the Riemann surface of y consists of a contour joining ot\ to 
x on the left of the cut in the first copy of C and a contour joining a± to x on the right of the cut on the 
other copy of C. Therefore 

G 1 (x)+ = G 2 (x)^ + f ydx . (4-13) 
From this the first assertion follows immediately with a = —i J ydx. Q.E.D. 
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Note that h(x) — diG 1 (x) and has no discontinuities (as we know already) on the gaps S J ; where it 
is actually harmonic, and is continuous on the cuts S/s where -however- is not differentiable (but it 
admits normal derivative on both sides). 

From the definition of Stokes' ray (4-3), we can define sectors Sk in U = C/(Uf +1 Ej) such that 

1. As x — > oo, the sector Sk is bounded by lZ 2 k and TZ 2 k+i (4-3) 

2. Let the y+ be the set 

y+ = {x € U, h{x) > 0} 

then each sector is asymptotically contained in y + . 

We will need later the following counting of edges and vertices 

Lemma 4.5 Let B = LlfL 1 Bi be the decomposition of the branchcut structure into connected components. 
Let bi be the number (which is even) of vertices and the number of edges in Bi. Then 

e; = 6,-1 (4-14) 

Therefore the total number of edges E is 

E = 2g + 2-K . (4-15) 

Proof. The first formula follows immediately by induction on the number of vertices; the second follows 
by considering that the sum Yli=i h — 2<? + 2. Q.E.D. 

4.2 The Stokes— Kirchoff normalized differential of second kind 

A crucial role will be played by a suitable normalized differential of the second kind which we baptize 
"Stokes" because of its relation (to be shown) with the Stokes matrices. 

Definition 4.5 A Stokes differential is a second-kind differential dw with the following properties 

1. it is antisymmetric w.r.t. the hyperelliptic involution; 

2. it has poles (without residues) of degree at most g + 1 at oo±; 

3. all periods around the connected components of the branchcut structure B are zero (or integer 
multiples of2iw). 

Some properties follow immediately: 

• The Stokes differential belongs to a Z K ^ 1 lattice of affine spaces of dimension 2g + 1 — K. Indeed 
there are a span of 2g antisymmetric second-kind differentials with poles of that order at oo± and 
without residue; imposing that the contour-integrals around each of the connected components is 
a multiple of 2in (the residue at oo± is zero so there are only as many loops as the finite gaps) 
imposes K — 1 affine constraints. 
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• Since the periods are zero (mod 2z7rZ) around Bi's, its Abelian integral w(x) = f£ dw defines a 
single valued (mod 2inZ) function on C \ B with a pole of order at most g + 1 at infinity. 

• By using similar arguments used for the G- function (together with the antisymmetry), we can show 
that it solves a certain scalar Riemann-Hilbert problem 

w{x)+ = — u>(x)_ + ifij , x 6 Sj (4-16) 
Mi -0 

on the edges of B. Here Si is the edge attached to a\ (the Weierstrass basepoint of the integration 
defining to); \i\ is zero because ot\ and x <G Si can be joined by a contractible loop whose projection 
on the x-plane intersects S only at x and ot\. 

The other important condition that we need to impose is specified in the next definition; note first 
that since there are no closed loops in the graph B we can always choose orientations of the edges of B 
in such a way that all the orientations of edges incident to any vertex are all incoming or outgoing. We 
assume one such orientation in the following definition. 

Definition 4.6 The Stokes differential dw is said to satisfy the (complex) Kirchoff's law if for each 
trivalent edge of B, denoted Si 1 ,Sj 2 ,Sj 3 the three incident edges (with their chosen orientation), then 

e i Wl + e t^ 2 + e » W3 = q ( 4 _ 17 ) 

where all edges are oriented so that at each trivalent vertex they are all incoming or outgoing. 

It should be clear that such Kirchoff-Stokes differentials do exist, since the Kirchoff's constraint poses 
only some number of nonlinear constraints. We need the exact count of these constraints. 

If hi is the number of vertices in Bi then there are y + 1 univalent vertices (of type I) and hence ^ — 1 
trivalent vertices (recall that hi is even). Summing up over all connected components of B we find the 
total number T of trivalent vertices 

t =e(|- i )=»+ 1 -^ ( 4 - 18 ) 

Therefore 

Proposition 4.1 The Kirchoff Stokes differentials on a Boutroux curve form a Z 2g+1 - K lattice of 
manifolds of dimension g = genus(C). 

The lattice aspect is due to the obvious fact that we can arbitrarily add integer multiples of 2tt to each 
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4.3 Asymptotics 

The following lemma shows that the spinors constructed in section 3.2 satisfy a Ricmann-Hilbcrt problem 
closely related to the one satisfied by G. This will allow us to express the large N asymptotics of certain 
orthogonal polynomials in terms of these spinors. 

The relation requires that we specify the choice of the flat line bundle Z v of Sec. 3.2 associated to a 
Stokes-Kirchoff differential: we will set 

rj = dw + Nydx . (4-19) 

Note that this n has residue ±N £ 1 at 00 ± (because of the normalized Boutroux condition). We will call 
the line bundle associated with this choice of Stokes-Kirchoff differential the asymptotic line bundle 
and denote it with 2,^. 

We note that it is the tensor product of a line bundle £,dw with the line-bundle 2,Ny&x] this last one, 
because of the Boutroux condition and since N £ N, is a unitary line-bundle. 

Lemma 4.6 Let £ be defined as in Sec. J^.l and the half integer characteristics chosen as specified in 
Prop. 3.1 for the given A and choice of orientations of edges of B 13 . Then the matrix 

h 
— ic r -i 



W := #jv, r = 



E r (x)e^ NG+w ^ as (4-20) 



(where the constants c r are specified in the proof: see eq. (3-68) for the definition ofa r (x)) satisfies the 
following Riemann-Hilbert problem 

( g-iWiTj-vj \ 

*+(*) = *-(a0 ( _ e iNa j+ i H o y ze£j 

*(a?) = of- l — r] , x^a 3 

tf(a;) = (l + C(x- 1 )a; rCT3 , x -> 00 (4-21) 
and ^(x) is holomorphic on C = C \ E. 

Remark 4.1 ^4n essentially identical RH problem appears in the work [12], where it was similarly reduced 
to the problem of Prop. 3.1. The idea perhaps was first used in [13]. 

Proof of Lemma 4.6. The growth condition near the branchpoints follows from (3-83). 

From the general definition in eq. (3-31) specialized to the differential (4-19) we see that the matrix 
^ is analytic in a punctured neighborhood of 00. Indeed the essential singularity of the spinors ir r , ir* at 



13 The orientation of the edge Ei attached to the basepoint of integration forces the orientation of the edges of the 
connected component of B to which it belongs (all edges should be either incoming or outgoing from the trivalent vertices). 
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oo is removed by the multiplication by the exponential exp J 77(73. This leaves us with a power growth at 
00; specifically 

nMx))^ Ny&X+dW 



'dx 



(4-22) 



has a pole of degree r at oo + (which is on the sheet chosen for p = X^ 1 (x)) and a zero of order r at oo_. 
This is so because (by their definition) the spinors ir r ef 11 have a pole of order r + 1 at oo + and a zero 
of order r — 1 at oo_; since dX has a double pole at both 00 ± and appears in the square-root, (4-22) 
follows. Similar considerations show that 

■ /^ (x) Nydx+dw 1 

(4-23) 



<(p(z))e 



x — > OO . 



Vda; c r 

Formulas (4-22, 4-23) define the constants c r ; the reason why c r appears in the denominator in (4-23) is 
simply that -by construction of Serre duality- 

1 



OO 4. 



res ( c r x r + 0(x r ~ i ) ) ( —x~ 
c r 



0(x- r ~ 2 ) )dx 



(4-24) 



If we remove the exponential part from W then it satisfies the Riemann-Hilbert problem (3-81, 3-82); 
with the exponential part the character specified by 

exp £ 77 : tti (C) - 

7 1 



§^ Nydx-\-dw 



(4-25) 



appears as the jump relations of the Abelian integral J Nydx + dw, thus yielding the proof. Q.E.D. 

For completeness we report the expressions for the entries of the matrix \&jv> ■ The two vector complex 
characteristics e N and (—S N ) are half of the a and 6-periods (on the chosen dissection of the curve in 
terms of a homology basis) of the differential Nydx + dw following the notation of eq. (3-31) 14 : 



Pa 



(Nydx + dw) 
(Nydx + dw) 



(4-26) 



After straightforward simplifications, the entries of ^N.r read (using (3-73) specified to the case of 
line bundle associated to 77 = Nydx + dw) 



-Q A (X(p)) , 
1 

-Qa(X( P )) 



„ e^( P -oo_) 


~A+e 
B+5 


N 

v _ 




- roo_ - 


-(r + l)oo+) 


^e^ 1 (?-«>+) 


e 


'A+e N ~ 


(r(oo- 


-oo+)) 



e r A ( P -oo + ) 


-A- 
—B- 


■*N. 




- roo + — (7* + l)oo_) 




e 


-A 
—B 


- S N. 


(r(oo+ — oo_)) 



*We only add a subscript N to emphasize the dependence on the large parameter N. 



36 





' A+e 
B+S 


N 

N _ 


(pH 


- (r — l)oo_ - 


roo + ) 


Cr-l Q r A {p-00 + ) 


e 


'A+e N ~ 


(r(oo_ — oo^ 


)) 



(*^.r) 2 i = ~— or .1 "/ i , ; Qa(^(p)) , 



1 e^- 1 (p-oo + ) 


-.4- 
— B- 


" e jv 


(pH 


- (r — l)oo+ - 


- roo_) 


Cr-l @a(P~°°-) 


e 


-A 
—B 




(r(oo_|_ — oo_ 


)) 



(*N,r) 22 = —^7^ T— r , . n . . " Qa(X( P )) , (4-27) 

f. r+1 

c r ■= r , (4-28) 

fc)^(oo + — oo_) 

where C\ was introduced in (2-18). The constants {c r ]>ez that appear in the formulae above have been 
determined by the normalization that ^fi t i = x r (l + 0(l/x)) and ^2,2 = x~ r (l + 0(l/x)). 

4.4 Semiclassical generalized orthogonal polynomials 

The main theorem of this section is that -broadly speaking- when (£, x,y) is an admissible triple the 
spinors ir n associated with it give the large N asymptotics of some orthogonal polynomials of the type 
considered in [2]. 

These polynomials are defined and related to our Boutroux-admissible curve as follows: let (£, x, y) be 
an admissible triple as in Def. 4.2 and let us choose a Kirchoff-Stokes differential dw (and corresponding 
Stokes function w) as in Def. 4.5 and with the same notations as in the remarks that follow that definition. 

Each component Bi with bi vertices has 6^/2 + 1 type-Ill (univalent) vertices on the boundary of as 
many distinct components 3^' of 3^+, each of which contains (at least) one Stokes sector; we can choose 
bi/2 + 1 distinct oriented contours T^e, £ — 1, . . . , y connecting one Stokes sector to the remaining ones 
and passing through the edges of Bi. It is easily seen that each edge belongs to at least one such contour. 

We associate a complex weight to each oriented contour I\^ as follows; the periods of dw which 
enter its RHP (4-16) define "complex currents" 

p 3 := e~ 1 ^ , (4-29) 

on each edge Sj of B. Because of the Kirchoff 's condition in Def. 4.6 we can always find other complex 
currents to associate to the contours Tij in such a way that the "net current" through each link 
is precisely pj. Note that p\ — 1 always 15 (sec for example Fig. 5): note that the weights n^i do not 
depend on the 27rZ arbitrariness entering the definition of the \Xj§. 

We repeat this procedure for all connected components B\ , . . . Bk of the branchcut structure B. A 
visual example is contained in Fig. 6 and its detailed caption. 

Define now the integral operator 

K bj2 

EE**.</ • ( 4 - 3 °) 

1 = 1 1 = 1 Jr i.' 



15 This will mean that one of the Stokes matrices is normalized, which we can always accomplish by a conjugation. 
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Figure 5: A connected component Bt with bt — 6 vertices, bi — 1 = 5 edges and bi/2 — 1 = 2 trivalent 
vertices. Shown are the net currents pi,.. . ,ps through the edges, satisfying the 2 Kirchoff constraints; 
we can find always appropriate complex currents to the contours (in textured linestyle) so as to give the 
desired net currents through the edges. 

Recall that the Boutroux curve defines a "potential" V(x) via (4-1). Let p n (x) be the monic orthogonal 
polynomials such that 

i I dx Pm (x) Pn (x)e~ NV ^ = h n S nm . (4-31) 



Denote by 4> n {x) and 4> n (x) the following functions 



±NV{x) r -±NV(z) 6 I \ 

Am J ^ z — x 

With these preparatory notations and remarks we are ready to formulate the main theorem of the 
section. 

Theorem 4.1 (Strong asymptotics) We have the following asymptotic estimates 
• For any compact set not contained in C\B, the following matrix-valued function 

* / x ( (t>N+r{x) (j>N+r{x) \ 

\ h N+r _ 1 < PN+r^l(x) Wl fa+r-l(l)y 

has the asymptotic behavior as N — > oo and r £ Z fixed 

t>N,r ~ y N Ax)e- NG{x)a3 , x e C/(Ufl 1 1 S J ) (4-32) 

away from the branch cuts. The matrix ^N,r was introduced in (4--20)) and the entries given by 
(4-27). 
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Figure 6: In this figure (which is only suggestive and not a numerical output) the grey area is where 
h < and the black thick lines are where h = and constitute the branchcut B; here the function h is 
continuous but not harmonic. The thin line in the grey area denotes the lcvelset of h(x) — c < passing 
through a saddle-point. Since there are 9 sectors within as many white areas the picture corresponds to 
some potential V(x) of degree 9. One should conveniently think of the grey area as the "sea" and the 
white areas as the emerged "continents". The paths Yk.j should connect the infinity within a continent 
to the infinity within another continent never going on the sea; for this purpose they can pass through 
the branchcuts as "causeways" . In each (oriented) path runs a complex current Xk j and the net currents 
within the edges of B automatically satisfy Kirchoff 's law. The only genericity requirement is that none 
of the net currents ought to vanish. In this picture there are only two components B\ and B 2 (left and 
right respectively) with 10 and 4 branchpoints; B\ connects 6 continents, each of which contains a sector 
where $l(V) — > +oo, whereas Bi connects 3 such continents. The paths T^j (dashed and dot-dashed) 
once defined can be deformed arbitrarily respecting the connectivity they provide between Stokes' sectors 
without affecting the values of (4-31) due to Cauchy's theorem (they are already depicted after the 
deformation but the reader should think of them at first as going exactly through the thick links). The 
existence of a situation with all the topological features of this example (and many others) will be proved 
in the last part of the paper. 
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On each edge E^ of B the asymptotic behavior as N — ► oo is given by 



/ / e -JVG(x) o \\ 

$N,r~ [^nA®) y ±p .-l e NG(x) e NG(x) j j > X £ E, (4-33) 

Proof of Thm. 4.1. From [2], the matrix Qn^(x) satisfies the following Ricmann-Hilbcrt problem 

(® N ,r) + (x) = (*N,r)-(.x) (l Ki f), r- I : |J I , 

(x) exv (^NV( X ) + (N + r)lo gX \ 

v N , r {x) - exp^ Q -(N + r) log a: J ' 

where the + and — indices denote the values of the function on the left and right hand sides of the 
contour respectively. The piecewise constant function k{x) is given by 

k{x) = { XeT * AB K (4-34) 

We can now transform the above Riemann-Hilbert problem by multiplying e NG ( x ) a 3 on r jght of 
(&n(x), where 03 is the Pauli matrix 03 = diag(l,— 1). Let 

&(x) $ w (x)e WG(K)ff3 . (4-35) 

This new matrix has jump discontinuities on E and T (see Figure 4 for E). The Riemann-Hilbert 
problem is transformed into the following for <5> 1 (x): 

§\(x) = $i(x)f v e V(G_- G+ ) ), x e r d s 

/ e ^(G+-G_) \ 

& + (x) = t>l(x)( e N(G-G + )), xeE/(snr) 



^(x) = [I + O 



x 



r<j 3 



Note that E \ (E fl T) is the union of all gaps and there is actually no jump on the last gap that extends 
to 00 (because the residue of rj at infinities is an integer) . 

Single cuts. Suppose Bi is a connected component of B with only one branch cut Ej (and hence 
Bi = £«) 

The jump matrix of $ (x) on the branch cut Ej is then 



*l, ^/e^ G +- G -) Pl e- N ^-+ G +y\ 



As in [14] and [33], we will make use of the factorization 



a -b \ _ / 1 W -b\ (I 

a- 1 -{ -a-'b- 1 l ib- 1 Oil -ar 1 1 
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and write the jump matrix as 



e JV(G+-G_) pie -N(G-+G+) 

gJV(G_-G+) 

1 0\ / Pi e- W ( G -+ G +) 

--IgiUVG- ! J I r l e iV(G_+G + ) o M r l e 2iVG + x 



(4-36) 
(4-37) 



We can now follow the same technique in [14], [33] and deform the contour into a 'lens' as follows. Let 
rf and be contours joining a^i-i and ot^i on the left and right hand sides of £, such that /i(x) is 
negative on both Yf and Yf (See figure 7) . This is possible because of admissibility. 

We can then deform the Riemann Hilbert problem near these branch cuts as follows. Let 

<& 2 (x) — <& (x), for x outside the lens-shaped regions, 

Q 2 (x) = ^(x) f p-i^2NG{x) l^j 7 for x in the left hand side lens regions (4-38) 
$ 2 (x) = <f> 1 (x) ( ^-i e 2JVG(x) j J ' f° r ^ m the right hand side lens regions 




if 

Figure 7: Lens-opening near an isolated branchcut. 



Multiple cuts. Now suppose that Bj is a component of B that contains more than one branch cut. 
Then apart from the boundary points of Bj, all the branch points are of type I. Let the boundary points 
be ae 1 , . . . , ai n . Since this connected component of B is a loop-free trivalent tree (in a finite region of 
the plane) and a neighborhood of it (less Bj) lies in the set h(x) < by definition of admissibility, we 
can then join these points by curves that lies within {h(x) < 0} as in figure 8. 

Let the curve between and dg i+1 be r^. From each trivalent vertex we run three contours joining 
the vertex to the closest amongst the arcs r^s without intersecting the cuts already made (see Fig. 9 for 
a self-explanatory exemplification) . On each side of each edge S m of Bj there are precisely two regions 
bounded by S m a contour amongst the r^'s and the added contours from the trivalent vertices; we denote 
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these two regions -D to .l and D m ^n; we then deform the Riemann-Hilbert problem near Bj to an equivalent 
one for a 2 x 2 matrix- valued function Q 2 (x) defined as follows in the various regions 

$ 2 (a;) = < & 1 (a;), for x outside the lens-shaped regions, 
>S> 2 {x) = $\x) ( __i 1 2 jvg(x) ?), xeD m<L 

$ 2 (x) = * x (a:) ( - lp 2 NG{x) ? ) , i£ £> m ,.R 

for each branch cut S m S i3j. 

The matrix <& 2 (x) then has jumps on F U X and the contours shown in Fig. 10. On (T U X) \ B it has 
the same jumps as < & 1 (x) (eq. 4-37). On the contours in Fig. 10 its jumps are the following: 

V 2 {x) = (-l e 2NG(x) j^j > x e 

= ( p -iJirG<*) J), (4-39) 
^ 2 0) = ( / r^NGM ? ) ! ^ e r mi (4-40) 



-iV(G_+G + ) 





where f m i(x) is either of the following, depending on the orientation of the contours 

f m i{x) = ±p^ l 1 ±pl 1 

in either of the above case, fjk(x)e 2NG ^ — > as N — > oo uniformly away from the branch points 



Branch cuts 




Added curve 



Figure 8: Lens-opening near a tree of B. 
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Branch cuts 



Added curve 



Figure 9: Lens-opening near a tree of B and added jumps. 




Figure 10: Another tree and lens-opening with added jumps. 
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From the jump conditions (4-39) and (4-10), we see that, away from the branch points, the Riemann- 
Hilbert problem for <I> 2 (x) approaches the following uniformly as N — > oo. 



i> 2 (x) $°°(x) 



where <&°°(x) satisfies the following jump conditions 




) 



—iNcTj 







x e 



) 



X G Sj 



(4-41) 



and <j>°°(x) — > / as x — > oo. 

By lemma 4.6, the matrix 'I'jv.r satisfies such the Riemann-Hilbert problem (4-21) and has the same 
behavior as x — > oo. Hence <&°°(x) = \&jv>(x). The asymptotic formula (4-32) and (4-33) can now be 
seen by reversing the sequence of transformations. To complete the proof, we need to show that there 
exists parametrices near the branch points. We will discuss this in the next section. 

4.5 Parametrix near branch points 

Since the Riemann-Hilbert problem for <i> 2 (x) does not tend to the one satisfied by < J>°° (x) near the branch 
points, we still need to solve the Riemann-Hilbert problem exactly near the branch points. The solution 
of the Riemann-Hilbert problem (4-39) for <i> 2 (x) near the branch points are the local parametrices and 
they can be constructed from the Airy function. This problem has been discussed in many places in 
the literature [11], [19], [14] and we shall not repeat the details of the analysis here. Instead, we would 
illustrate the main idea and provide the solution to such a problem. The only relatively new feature in 
this respect- is the presence of turning point (the branch-points) of type I, namely where three branchcuts 
are connected. 

In order to complete the asymptotic analysis and obtain a uniform estimate, we need to find a solution 
<& p (x) of the Riemann-Hilbert problem (4-39) in a neighborhood O a of a branch point a such that as 
N — > oo, we have <& p (x) — > <& 2 (x) near the boundary of such neighborhood. We transform the RH 
problem a first time by multiplying the solution $ 2 (x) by exp(— NG{x)a^) on the right 



As a consequence, the Riemann-Hilbert problem satisfied by the matrix Z v {x) is transformed into a 
problem with constant jumps. 



Z p (x) = $ 2 (x) exp(-iVG(x)cr 3 ), x e O, 
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Univalent turning point (type III). If a is a turning point with only one incident branch-cut, say 
S m , then the Ricmann-Hilbert problem for Z v in O a consists of the following jumps on the contours 7* 



where p m is a constant depending on which branch cut is incident to the turning point under scrutiny. 
(See figure 11) 

Trivalent turning point (type I). If a has 3 branch cuts attached, the function Z p (x) then satisfies 
the following Riemann-Hilbcrt problem 



where the orientations of the rays ji are chosen such that all of them are pointing away from the branch 
point. 

A different choice of orientations will only result in a change of sign of some of these constants. The 
crucial observation is that with this choice of orientations we have p mi + p rri2 + p m3 = 0, which implies 
that the product of all the above jump-matrices is the identity and is precisely the Kirchoff 's law for our 
"currents" pj's. 

We would like to construct solutions Z p (x) to the Riemann-Hilbert problems (4-42) and (4-43) near 
the branch point a such that 




(4-42) 




(4-43) 



Z p (x)e 



,NG(x)a 3 



$°°(x)(l + 0(N- 1 )), xedO a , N^oo 



(4-44) 
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The solutions near these branch points can be constructed by using the Airy function. The Airy function 
Ai(x) is the unique solution of the ODE 

with the following asymptotics 

M() . _l. r ^( 1 + (_L)) 

M® - -^ { .«-Kl( 1 + (^)) 

as £ — > oo and | arg£| < ir. The branches of £J and in the above are principal branches with branch 
cut at the negative real axis. 

Since G(x) — G(a) behaves like cx i near a branch point for some constant c, we can choose a function 
f = (%N(G(x)-G(a)))s that maps the branch cut onto the negative real axis. This is possible because 
the function %N(G(x) - G(a)) maps the branch cut onto the imaginary axis. The function £ then defines 
a one-to-one mapping between the neighborhood O a of the branch point a and a neighborhood of the 
origin in the complex £ plane. (See Figures 11, and 12.) 

7i 

7i O a \ 





Figure 11: The function (|iV(G(x) — G(a)))* maps the neighborhood O a of a to a neighborhood of the 
origin. 



As in [11], a solution of the Riemann-Hilbert problem (4-42) is given by the following. Let ^f p (x) be 
the matrix given by 

ZP(n _ (MO M^ 2 xeI 

16 This requirement is due to the fact that in defining the Airy function Ai(£), we assume | arg£| < it. The requirement 
can be dropped if we use different branches of the Airy function. See [19] 



4G 



o a 




-73 



7i 



Figure 12: The function (|iV(G(x) — G(a)))2 maps the neighborhood O a of a to a neighborhood of the 
origin. 



z p (0 = 
z p (0 = 
z p (0 = 




(4-45) 



where u> — e 3 . Then by using the relations 

Ai(£) + uAi(w£) + uo 2 Ai{& 2 ) = 
Ai'{£) + Lu 2 Ai'(£uj) + uAi'{& 2 ) = 

one can show that Z p {£) defined by (4-45) does satisfy the jump conditions (4-42). 

Since the only jump discontinuity of £5 is at the branch cut where it changes sign, while N(G(x) — 
G(a)) has jumps at the gap S TO and the branch cut E m , these two functions are not identical. However, a 
branch of £2 can be chosen such that the difference between |£s and N(G(x) — G(a)) remains bounded 
in O a . 

Let £ be the function NG(x) — where the branch of £5 is chosen such that £ remains bounded 
in O a : in fact this difference is a locally constant function on O a \ E. 

By considering the jump discontinuities of NG{x) and ^ , we see that £(#) is bounded throughout 
O a and has the following jump discontinuities 

((%)+ = C(x)- + iNa m , x € S m n O a 

C(x)+ = -({x)_+tNa m , x€£ m nO Q (4-46) 
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To fix the boundary conditions such that Z p (£ i )e NG ^ c ' 3 — > $°°(x) as — * oo on the boundary of O a , 
first note that Z p (^)e NG ^ (T3 has the following asymptotic form 

Let Z a s(<£) be the leading term of the above asymptotic expansion. Then from (4-46) and the fact that 
£i — > i£i across the branch cut of G(x), we see that Z as (£) has the following jump discontinuities 

Zas{£)+ = Zas(£)- ( _p-l e iNa m q )) £ € £ TO PI O a 

/ e iNa m q \ 



Therefore the matrix 



i^x)^ 00 ^" 1 ^) 



is single- valued in O a . It is also holomorphic in O a because it can at worst have a square-root singularity 
at a, but since it is single- valued, this cannot happen and hence E(x) is holomorphic in O a . 

Hence the function E(x)Z p {x) is bounded in O a , tends uniformly to $°°(a:) near the boundary of O a 
and satisfies the jump condition (4-42) . 

For a branch point with 3 branch cuts attaching to it we choose £ = (%N(G(x) — G(a)))s that maps 
the branch cut £ mi onto the negative real axis. The following Z p (£) would satisfy the jump conditions 
(4-43) 

»® - {% JP&y^ ?)Gfifc)'. ««' 

- -(;?</) ^)J-"l : . iJUSfc)*' 

- -(-1 -mS?)— (- , J)te)* {6 " <" 7 > 

By using p mi + p m2 + p ni3 = 0, one can check that Z p (£) does indeed satisfy the jump conditions (4-43). 
Let be the following function 

C(x) = NG{x)-£, x e IUIIUVUVI 
C(x) = iVG(x)+£5, xelllUlV 



Ai'(0 lo 2 M'{^OJ \p miPm2 
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where the branch of £2 is chosen such that C{ x ) is finite in O a . Since £2 has only one branch cut at 76 
while NG(x) has 3 branch cuts at 72, 74 and 76 respectively, £(2) has the following jump discontinuities 

C(x)+ = -C(x)- + iNa mi , x e 76 

C(x)+ = -C(a:)- + iNa m2 , xe l2 (4-48) 
C(x) + = -C(z)- + iNo m3 , x € 74 

To fix the boundary condition in this case, first observe that the function Z p (x) has the following asymp- 
totic expansions 



V —1 1 / V V N 



CT 3 



2 V 7r \ _1 V \ \ N J J \P mi Prn 3 

as £ — > 00. If we now take take Z a s(£) to be 

ITT ^13. 

2\/7r V -1 V \p mi Pm 2 J 



ITT 53 

2\/7T V -1 1 / \PmtPm3 J 



then from (4-48), we see that Z as (x) has the same jump discontinuities as <& 2 (a;). 
Therefore the matrix 

E(x)=^(x)Z-J-(x) 

is holomorphic in O a and E(x)Z p (x) will be bounded in O a and will tend uniformly (to within 0(1/N)) 
to <&°°(x) near the boundary of O a . Also, it satisfies the jump condition (4-43). 
Finally, if we let <& 3 (a;) be the following matrix- valued function 

$ 3 (a;) = $ p (x), zeug! 2 Q , 

where $ p (a;) = E(x)Z p (x)e NG ^" 3 . Then as in [11], [14] and [19], one can show that the Riemann- 
Hilbert problem satisfied by Q 2 (x) tends to the one satisfied by <J> 3 (cc) uniformly at all points as N — > 00. 
Therefore the function <& 3 (x) would give the asymptotics of the orthogonal polynomials as N — > 00. 
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4.6 Density of zeroes of the orthogonal polynomials 



The content of Thm. 4.1 is that the (monic) OPs pN+r behave uniformly on compact sets not intersecting 
the branchcut structure B as the following expressions (see 4-27) 



p N +r(x)e ™ v{x) ~ (*JV,r)l,ie 



Nq(x) ~ @a(P-°°-) 6 


' A+c 
B+S 


N 
N _ 




- roo_ — (r + l)oo+) 


r e^ +1 ( P -^ + ) 


e 


'A+e N - 


(r(oo_ — oo+)) 



Q A (X(p))e-^^ 



where c r is the constant (independent of N) computed after (4-27), p is point on the curve C with 
X(p) = x and on the same sheet as oo + and the vector of complex characteristics e N ,6 N are given in eq. 
(4-26). 

In particular (using that the Abel map of oo + is opposite of that of oo_) it follows that 



N 



In \p N+r (x)\ 



1 

2 

1 

N 



$t{V{x)) - h(x)- 
ln 



Q A (x(p))e r A ( P -^) e 


1 S n\ 


(ph 


- rcx)_ — (r + l)oo_|_) 


e^ip-^) 


e 




(r(oo_ — oo + )) 



In \c r \ 
N 



, /i(x) := 5R 5 (x) 



(4-49) 
(4-50) 



uniformly over compact sets in C \ B. 

In the RHS of (3-73) there are at most g zeroes (on the first sheet) which do not belong to B; therefore 
we have immediately 

Proposition 4.2 All but at most g zeroes of the polynomials pn+i- are contained in an arbitrary neigh- 
borhood of B for N large enough. 

In different terms this means that the almost all zeroes accumulate on B. 

It should be evident that the zeroes do not actually lie on B but they get closer as iV increases. From 
Thm. 4.1, formula (4-33) it follows that in a neighborhood of B we have 



PN+r (x)e~ NV ^ ~ Ae NG ^ - Be~ NG ^ , 



(4-51) 



where A, B are expressions depending in x but with at most g zeroes (whose count is asymptotically 
irrelevant) . 

Solutions are of the form 



±^^{B/A) + 2i^ 



(4-52) 



which shows (in a slightly heuristic way) that the zeroes have asymptotically vanishing -KG (which means 
that they are "close" to the branchcuts B where h = KG vanishes) and A3G ~ -k ■ 
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Using Cauchy-Riemann's equations for G one finds that the tangential density along a cut is 



i d$G i dh 

— «— = 75 — 4 " 53 

IT OS 7T on + 

where ds and dn denote the tangential and normal derivatives along the (smooth parts of the) branchcuts 
B. 

To put it differently the density of zeroes of pn+t per unit length along an arc of B as N — > oo tends 

to 

*-M-s(£ + si)~i£- » S \UW (««> 

where <9n-|- denote the normal derivatives on the two sides of the cut. 

Note that this expression is positive because the function h is zero on the cut and negative (strictly) 
on its left/right neighborhoods (by definition of admissibility of our triple (C,x,y)), hence the normal 
derivatives are strictly negative on the cuts (but away from the turning points) . By elementary harmonic- 
function theory (or using electrostatic analogy) we then have 

h(z) = InV(z) + f In \z - C|^oo(C)d S (4-55) 

z JB 

(where ds denotes the ordinary arc-length) and hence the total mass of J„ u 00 {x)ds is necessarily 1 since 
that corresponds to the residue of ydx (and to the constant in front of the logarithmic term of G(x)). 

5 Part III: Reconstruction of (admissible) Boutroux curves 

In this section we prove that any admissible graph is the graph of an admissible triple. In order to do 
this we take a detour in the theory of Strebel differentials, which we state here in a simple form suitable 
for our application; all the general statements can be found in [32] and [22]. 

Let P(i 
differential 



Let P(z) — n!Li( z ~ a i) tI ' i be an arbitrary polynomial of degree n = It defines a quadratic 



<P(z) := P(z)dz 2 . (5-1) 

The metric associated to a quadratic differential <fi, denoted by |</>|, is (quite literally) |P(z)||<iz| 2 ; in our 
case this reads 

|^| = |P(^)|fd a; 2 +dy 2 V (5-2) 



It is an easy exercise to see that this metric is flat, with conical singularities at the zeroes of P(z) 
(and cuspidal singularities at the poles of order higher than 1). The flat coordinates arc given by the 
real/imaginary parts of 

w := J v/P(s)ds (5-3) 
which is a locally defined parameter away from the zeroes (poles) of P. 
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Definition 5.1 The horizontal lines are the lines defined by —const ( or equivalently arg(-\/ 'P(z)dz) = 

0,7T). 

Since \/P is denned up to a sign, so is w (and up to translations as well), but the notion of horizontality 
is well-defined. 

Definition 5.2 The critical horizontal lines are the horizontal lines which contain any of the critical 
points (zeroes) of P. 

The critical horizontal graph (critical graph for short) is the union of all critical horizontal lines. 

It is a simple check in a local coordinate that from a zero z — a of multiplicity [i there originate fi + 2 
critical horizontal lines, with relative angles of 27r/(/i + 2) at a. 

Clearly the critical graph of P consists of the union of a finite number of Jordan arcs. 

In the general theory of trajectories of quadratic differentials one may encounter trajectories that fill 
domains with non-empty interior, called recurrent trajectories: by definition a trajectory 7 : R — > C (C 
an arbitrary Riemann surface) is recurrent if it belongs to its limiting set 

L -r ■= f| UW*)} ■ ( 5 "4) 

tglsX 

Visually these are curves filling "ergodically" some region. This is not the case in the situation at hands 
because of the following 

Proposition 5.1 (Thm 15.2 [32]) No trajectory ray of a holomorphic quadratic differential on a do- 
main of connectivity < 3 is recurrent. 

In our case the quadratic differential is holomorphic in CP 1 \ {00} = C: therefore (as follows also from 
[22]), our critical trajectories can only connect two critical points or one critical point to 00. 
We first dwell a bit on the topology of the critical graph, X. 

The statements summarized in the following lemmas can be found (with different notation) in [22] 
but are not difficult to show directly; here we restrict the attention to polynomial quadratic differentials, 
but more general statements can be found ibidem. 

Lemma 5.1 The graph X contains no loop in the finite part of the plane. 

Let C \ X — UTj be the decomposition in connected components. 
Lemma 5.2 Each connected component Tj is simply connected. 

Lemma 5.3 Each (simply) connected component Tj has at most two boundaries, each consisting of an 
infinite piecewise Jordan curve. Each of the boundary components contain at least one critical point (since 
it is constituted of critical lines). 
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In each Tj we can choose one boundary critical point zq as basepoint for integration and define wj to 
have ^s(u>j) > 0, uij(z a ) = 0. 

Lemma 5.4 The function Wj — \J P(s)ds is a uniformization function from Tj to an infinite hori- 
zontal strip or to the upper half-plane. 

Lemma 5.4 shows that C \ X is the union of halfplanes and strips (topologically) ; the complex z- 
plane itself, C 2 is then realized as a union of closed half-w-planes or horizontal w-strips with appropriate 
identification of the boundary points; each of these halfplanes/strips has at least one marked point on 
each boundary. 

Definition 5.3 A marked half plane/ strip is a copy of > (0 < 3u> < i, respectively) together 

with a collection of marked points on boundary (at least one for each boundary) and up to horizontal 
translations. 

From this definition it is clear that a marked halfplane with K + 1 marked points (K > 0) has K real 
moduli, whereas a marked strip with K + 2 points (at least one for each boundary) has K + 2 real moduli 
or, better, one complex modulus and K real moduli. These are simply the differences of two chosen 
critical points on different boundaries (which gives a complex parameter with nonzero imaginary part) 
and the K remaining relative positions of the other marked points on the two boundaries (which are 
real) . 

Near z = oo which is a pole of order n + 4 for the quadratic differential 4>, the general theory [32] 
shows that 

1. any horizontal line approaches oo asymptotic to n + 2 directions forming relative angles 

2. any noncritical horizontal line in a neighborhood of the pole is a topological circle, approaching the 
pole along two consecutive asymptotic rays. 

This in particular forces the following 

Lemma 5.5 The left of the two rims of any strip approach 
oo asymptotically along the same direction. Ditto for the 
right. 

Proof. Suppose that -say- the right upper/lower rims ap- 
proach infinity along different critical directions; then, by 
2 above, there would be a noncritical horizontal trajectory 
within the strip that is topologically a circle with one point 
at oo and confined in 5R(w) << or di(w) >> 0. Such a 
trajectory would not span the whole strip from one side to 
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Figure 13: The structure of the horizontal 
lines near infinity (in the local parameter 
( = -), here for a polynomial of degree 2. 




the other. This is a contradiction, since horizontal noncriti- 
cal trajectories within a strip span the whole strip from left 
to right. Q.E.D. 

Lemma 5.6 Two strips cannot have a complete boundary in 
common. 

Proof. If this were the case there would have to be at least one critical point on the separating boundary; 
such point would have only two critical lines, but then it would be necessarily regular (the critical points 
are all at least trivalent). Q.E.D. 

We can depict the topological structure of the graph X as follows: we draw a disk, whose boundary 
represents the asymptotic directions at infinity. On the boundary we mark n + 2 points representing 
the asymptotic directions of approach. Inside the disk we mark the critical points and connect them 
according to the connectivity dictated by the graph X (see the example in Fig. 14). 

It follows from Lemma 5.5 that the two left (right) rims of a strip approach oo along the same direction, 
therefore they join at the same vertex on the clock diagram. 

Decorations. Each strip has at least one critical point on each rim; we chose two such points a, (3 and 
associate to this pair the integral p :— y/P(x)dx, where the path of integration lies within said strip 
and the branch of the square-root is chosen so as to have a result with positive imaginary part. We call 
this number the modulus of the strip (associated to the given choice of pair of points). 

If the boundary between regions has more than one critical point belonging to it then we assign to 
each arc between two adjacent critical points the Strebel length of that arc, namely \J P(x)dx, where 
the integral is performed along said arc and the branch is chosen so as to have a positive result. 

The clock-diagram, together with these parameters (the complex moduli of the strips and the positive 
Strebel lengths of consecutive critical points on the same boundary) will be called the decorated clock 
diagram. 

The critical graph X can be considered (from a topological point of view) as a loop-free forest made of 
vertices of different valencies and edges connecting vertices with either other vertices or infinity (along a 
given direction). In the following lemma we study some elementary cnumerative properties of this graph. 

Lemma 5.7 Let /Ltj + 2 be the valencies of the vertices, V the number of vertices, P the number of 
halfplanes and S the number of strips. Let oi be the number of open critical lines from the i—th vertex 
(i.e. lines that go to infinity), and Ci the number of critical lines that are closed (i.e. go to some other 
vertex). Then 




(5-5) 




(5-6) 



As a consequence, the number halfplanes is n + 2 = 2 + J^Mi- 
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Figure 14: An example of a critical graph for a polynomial of degree 4 and its "clock" diagram. 



Proof. First of all Oj + Cj = fa + 2 because this counts the valency of each vertex. The sum ^ Cj is even 
because each closed line appears in exactly two vertices (i.e. is counted twice). Note that each halfplanc 
has at least one critical point (vertex) on its boundary, and there must be two open critical half-lines on 
its boundary. On the other hand each strip has 4 open trajectories (two for each side) . So AS + 2P is 
the number of sides of open trajectories, i.e. twice the number of open trajectories 2 o,;. From this we 
have the first equation (5-5). 

The second formula is proved as follows. Assume first that there are no closed trajectories; in this 
case each vertex is on the boundary of some strip, unless there is only one vertex (in which case the 
formula holds trivially). This is so because the i-th vertex splits the plane into Vi sectors, one of which 
must contain another vertex. Then the two rays bounding this sector are one side of a strip. Each vertex 
then contributes one strip and hence the number of strips is V — 1. 

Suppose now that the i-th vertex has Cj closed trajectories; this means that there are Ci other vertices 
connected to this one (necessarily distinct, since there are no closed loops): they do not contribute per 
se to the number of strips. The formula follows. 

Finally, by substituting (5-6) into (5-5), one has P = 2 + fa = n + 2. Q.E.D. 

Collecting Lemmas 5.2,5.3,5.4, 5.7,5.6 together with the idea of the clock diagram we see that any 
clock-diagram must follow the rules formalized in the following definition. 

Definition 5.4 A (topological) clock diagram is a regular n + 2-gon called the clock (whose vertices are 
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"infinite") containing V vertices (the "finite vertices") of valencies [ii + 2 such that 

v 

8=1 

Moreover 

1. each edge connects two finite vertices or one finite vertex and an infinite one; 

2. there are no loops in the interior of the clock; 

3. each component Tj in which the interior of the clock is separated by the edges is topologically a 
marked halfplane or strip. 

4- each half-plane borders exactly one of the sides of the clock. 
5.1 Inverse problem 

We want to show that for any graph X giving rise to a decorated clock-diagram with the topological 
properties of Def. 5.4, there is a polynomial whose critical graph corresponds to the given graph. This is 
essentially an (interesting) exercise in welding [32]. 

In fact the problem is essentially solved and in the more general setting of rational quadratic differ- 
entials in [22], so we basically rephrase the contents of [22] adding the details that are needed for our 
application. 

Suppose we have a clock-diagram; by Def. 5.4 the interior of the clock is partitioned in regions Tj 
each of which is either a halfplane or a strip (topologically) . This implies that the given diagram tells 
us how to glue (topologically for the time being) half-planes and strips so as to have a simply connected 
topological space. 

We now choose an arbitrary decoration of the clock-diagram: for each component of a boundary with 
K + 1 marked points we assign (arbitrarily) K real positive numbers, representing the relative distances 
between consecutive marked points; for each strip, we choose two marked points on the two distinct 
boundaries and assign an arbitrary complex parameter p, 9(p) > to this pair. 

By virtue of this construction we will have decorated marked halfplanes/strips in such a way that 
the decorations on the boundaries (i.e. the relative distances of the marked points) match between 
neighboring Tj-'s. 

We think of these marked abstract halfplanes/strips 1^ as realized in copies of the iu-upper-half-planc 
and introduce a flat coordinate Wj for each halfplane/strip, normalized so that Wj vanishes at one of the 
marked points on one boundary. 

We form the topological surface X = LiTj/ ~, where the equivalence relation is the metric identifi- 
cation of the boundaries according to the coordinates Wj and the topological structure dictated by the 
clock diagram. 



5G 



5.1.1 Conformal structure 

The construction of the conformal structure follows [22], with some minor deviations. 

Our topological surface X is connected and simply connected since it is a model of the interior of 
the clock (and hence of the plane): the conformal structure is defined in an interior point of Tj by 
the coordinate wj itself. A neighborhood O of a critical vertex a of multiplicity fi intersects fi + 2 
halfplanes/strips which we simply denote by Li, . . . , L M+ 2- The local coordinate z in this neighborhood 
O is then defined by 

z{p) = e^wf^, peOnr<, w e = V / P(s)ds, £ = l,...,n + 2. (5-7) 

J a 

Near a smooth common boundary points between -say- Li,L2 of coordinates s, H, the coordinates are 
(w% — s) and ±(w 2 — s) on the two sides (respectively) , where the sign depends on the relative orientations 
of the strips/halfplanes at that point. 

It follows from these definitions of the conformal structure that (dwj) 2 lifts to a well-defined holo- 
morphic quadratic differential on the (now) Riemann surface X; this differential has zeroes precisely at 
the marked points and are of multiplicities fi. 

In order to conclude that it defines a polynomial, we need to compactify the surface and show that 
the quadratic differential has a pole at the point of compactification. 

Compactiflcation We first topologically compactify X using one-point compactification ( Alexandrov) . 
The ensuing topological space is connected, compact and simply connected. 

Let n = fa be the sum of all multiplicities of zeroes of the quadratic differential dw 2 which -as 
discussed above- is globally defined on X \ {oo}. 

We need to define a local coordinate £ at the compactification point oo; since a neighborhood of oo 
intersects all regions Tj's, we need to express the to-be-defined coordinate £ in terms of the coordinates 
Wj naturally defined in each region; moreover we must do so in such a way that boundaries of adjacent 
regions are mapped to the same line in the £-plane near £ = C(°°) = 0. 

Let Uoo be a neighborhood of oo and let Vj be the intersections of Uoo with all the domains L/s. 
To have a pictorial idea we can turn the clock-diagram inside-out so that the inside of the (n + 2)- 
polygon represents the point oo and the lines from the vertices the critical trajectories (see figure 15 for 
an example). 

Let us proceed clockwise on such an inverted-clock diagram; because of our convention ^s(wj)'s the 
orientation of two consecutive half-planes (the regions that border one edge of the inverted-clock diagram) 
must be opposite, independently on the number and orientations of the subregions of strips that may be 
incident on the vertex between the two planes. 
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Figure 15: An example of the process of compactification and the corresponding inverted 
at the point of compactification. 
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Starting from region V\ we can continue the differential dw — dwi to the region V2, where dw — ±du>2- 
If V2 is one end of a strip then we may have one or the other sign, depending on our choice of uniformizer; 
if V2 is another half-plane then necessarily we have a minus sign according to our observation above. 

Recall that there are n + 2 = 2 + ^ fij half-planes (edges) ; therefore after going around 00 once we 
will have monodromy sign of +1 if n is even, —1 if n is odd. 

In other words we have 

Lemma 5.8 If n = 2k then the differential dw is well-defined (up to overall sign) in Uqo \ {00}. 
If n = 2k + 1 then the differential dw is well-defined on the double-cover of Uoo \ {00}. 

In the case of even n (the case of interest to us) we can define a function w by integrating (one branch 
of) dw; this function is not single- valued in general but has additive monodromy 

w 1— > w + <j> dw = w + 2inf3 . (5-8) 

In fact it is not difficult to see that (3 has a geometrical meaning of the alternating sum over all vertices 
of the clock of the total modular width of the strips incident to each vertex. 

If n is odd it is also not difficult to see that § dw (where the integral winds twice around the com- 
pactification point) is always zero due to the monodromy dw 1— ► — dw after one loop. 

A uniformizer £ — C( w ) must have the following properties 

• C( w ) is a locally analytic function for \w\ >> 0; 

. c(^ + /?) = CH; 

• ((w) has a singularity of type ui~5 _1 (so as to be able to accommodate n + 2 halfplanes in one). 
We sketch the main steps of the construction following the ideas in [23, 32]. 

Consider the intersection of a halfplane, denoted by III, with a neighborhood of infinity; III comes 
equipped by definition with a coordinate w\ . We now proceed clockwise around infinity and consider the 
next region of intersection; this might be a strip or another halfplane which borders III on the boundary 
%{wi) =0, 5R(wi) » 1. 

Using now the coordinate of this region (possibly up to translations) we can extend the coordinate wi 
to this strip; if the other edge of this strip borders another strip we repeat the procedure until we reach 
a strip which borders with the next clockwise halfplane n 2 . 

The region we obtain looks like the one in Fig. 16 

Let's denote H Pl this region, where p± is the sum of all complex moduli of the parts of strips that we 
have attached to the right of the halfplane: we repeat this construction for all n+2 half-planes, obtaining 
thus regions H P2 , . . . , H Pn+2 ■ 
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w;-plane 



Figure 16: The regions used for compactification; on the left an example of a shifted region H p + p' and 
on the right the image through ((w) of that region; the asymptotic directions have an angular separation 
of 2Tr/(n + 2). Subsequent regions are mapped to the £ plane by choosing different branches so as to fill 
all sectors of a neighborhood of £ = 0. Also shown are the separations between the added half-strips on 
the bottom right of the w halfplane (in this case two strips). 

In order to glue one such region to the next, keeping in mind that consecutive half-planes have upside- 
down matching, we have to glue H pi with — H P2 + pi and so on and so forth. We obtain a chain of regions 
(we assume for definiteness n even) 

n+l 

H Pl -» -H P2 + Pl -> H Pa - pi + p 2 -> -H Pi + p x - pi + p 3 ■ ■ ■ H Pn+2 + ^(-l>Vj • (5-9) 

j=i 

The resulting Riemann surface V has in particular two boundaries corresponding to a half-line in the 
negative real w axis of H Pl and and a half-line in the line {w + p n +2, $s(w) = 0, $t(w) » 0} of the last 
region; we denote them by respectively. 

In addition on V the function w is now single valued and realizes an/2 + 1 cover of a neighborhood 
of oo; the two boundaries Ej, and S« are mapped to 

w(T, L ) c {3?H = 5R(w) « 0} (5-10) 

iw(Sfl) c {3(w) = 3(/3) 3?H « 0} (5-11) 

where 

1 n+2 

^=^E(-D^ (5-12) 

i 

and the points on these lines are identified by w ~ w + /3. Now consider the map ((w) defined implicitly 
(up to overall additive constant) by 

C"^ +/3ln( = w (5-13) 



60 



For \w\ large there is a univalent branch of C( w ) which behaves as any given branch of w~ n + 2 . 

Under this map, the line £l is mapped by (one branch of) C,{w) to a Jordan arc C(£l) approaching 
C = 0. Consider a neighborhood U of ( — with a cut along this arc, so as to have a simply connected 
domain. The image w(U \ £(Xl)) covers a punctured neighborhood of w — oo precisely n/2 + 1 times 
and realizes a biholomorphic equivalence with the Riemann surface V; because of the periodicity of Ci w )-> 
the two sides of the cut C(^l) are mapped to Xl and £_r. 

Thus the function £ realizes the above identification of and and gives a conformal uniformization 
of X; in addition £ is continuous at the compactification point oo and therefore defines a conformal 
structure at this point 17 . 

This construction also shows that the differential dw 2 has a pole of order n + 4 at the compactification 
point oo in the local uniformizer £. 

Since X is holomorphically equivalent to the Riemann-sphere the differential dw 2 is represented in a 
global uniformizing coordinate z as 

dw 2 = P(z)dz 2 , (5-14) 

with P(z) a polynomial of degree n: we have thus succeeded in proving the equivalence of polynomials 
and decorated admissible graphs. 

Remark 5.1 In a nutshell the above discussion boils down to the following statement: the space of 
polynomials of given multiplicities of zeroes 

K 

P(z) = Y[(z- aj r> , (5-15) 

can be (locally) parametrized by decorated admissible Strebel graphs. The "coordinates" are (essentially) 

Ej-i := f 3 ^P{z)dz , 3 = 2,..., (5-16) 

J ai 

up to translations and dilations that leave these integrals invariant (and for some choice of the contour 
of integrations and branch of the square-root). 

A naive parameter counting confirms this fact: polynomials of given multiplicities of zeroes are 
parametrized by the K positions of the zeroes (distinct) and an overall multiplicative constant, in to- 
tal K + 1 complex parameters. If we factor out the group of invariance of the above integrals, namely the 
action of translations z i— ► z + c and dilations 

z = \z, P(z) = \ 2 P(\z) (5-17) 

we have the same number of parameters, K — 1. 

17 If n is odd we have already shown that (on the double cover) w has no monodromy and the construction works identically 
with the caveat that we need to use the double cover of the neighborhood of the compactification point. 
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For our application, however, it is essential to control the global topological structure of the critical 
graph, something the above coordinates tell nothing about. One obvious reason is that the zeroes of the 
polynomials are on the same footing and hence -even for the same polynomial- we could assign different 
coordinates. More importantly, we could have two inequivalent clock diagrams with decorations given by 
the same numbers: the corresponding polynomials would then be different but with the same Strebel lengths 
between zeroes 18 . 



5.2 (Admissible) Boutroux curves 

In the case of relevance to our paper we must restrict the parametrization given by the decorated clock- 
diagram of a polynomial to a suitable submanifold; this is the submanifold of polynomials of the form 



2g+2 

P(z) = -M 2 (z) X[(z-a 3 ) = 



-y 



(5-18) 



(where the sign is just conventional for our application so that the function h introduced in Def. 4.2 is 
the imaginary part of the Strebel flat coordinate w) with additional constraints on the clock-diagram 
and on its decoration described below. 

The corresponding clock-diagram has only trivalent ver- 
tices (the a/s) and even-valent ones (the zeroes of M). The 
Boutroux condition (Def. 4.1), as discussed implies that the 



critical 0-level set 19 Xo of 3 J* ^P{z)dz is well defined 
independently of the choice of critical point otj and indepen- 
dently of the choice of contour of integration. 

In this case the height function h = 5w can be defined 
as a continuous global function with smooth saddle points 
at the zeroes of M(z): indeed it is easily seen that there 
is a consistent choice of signs for the Strebel coordinate near 
each even critical point (the zeroes of M all appear with even 
multiplicity in the quadratic differential) in such a way that 
the height function h = 3(ui) is harmonic there. 

The critical horizontal graph of P(z) thus naturally splits 
into Xo U Xi, where all the Qfj's belong to Xo- The distances 
between critical points aj belonging to the same connected 

18 To put it differently, clock diagrams with given number of finite vertices represent cells in the space of polynomials; 
the decorations give coordinates to each cell, but it is meaningless on a global level to consider only the decoration. 

19 The critical level set of a function is the union of all connected components of the level-set that contain at least one 
critical point. 




Figure 17: An example of reconstruction 
of admissible triple. The part within the 
graduated square is an actual numerical 
output . 
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component of Xo are arbitrary real numbers. The other mod- 
uli are the complex moduli of the strips: they must satisfy the real constraint that the heights of the 
strips "between" two connected components of Xo must add up to zero. 

The reconstruction theorems proved in the previous sections could be rephrased by saying "if you 
can draw the graph you want, then it exists (with the same topological features)"; instead of giving an 
abstract account of the procedure we find it more instructive and transparent to analyze some examples. 

Let us consider the examples in figures (17, 18, 20, 19). In Fig. 17 the black contours make up Xo, 
whereas the complete Strebel critical graph is the collection of solid and dashed curves. The numbers 
£i, £2 are arbitrary real (positive) numbers; they are f yPdx (where a, a' denote the two critical points. 
The two numbers pi are instead complex (with nonzero imaginary part) and correspond to f c yPdx, 
where c is the critical point in the middle. They satisfy the only constraint that 5(pi) = $s(p2). Of course 

there could be more strips between the two connected components of Xo- (The normalization res ydx — 1 

00 

imposes the constraint l\ + £2 = 1). 

The topology of this example is the "usual" one in the two-cut potentials: there are two arcs supporting 
the asymptotic distributions of zeroes of the orthogonal polynomials and the differential ydx has a critical 
point between them. The degree of the potential V is here quartic. 

The second example (Fig. 18) is "unusual" but it can happen for this class of semiclassical orthogonal 
polynomials: it is a cubic potential without double zero for the differential ydx (similarly, one could have 
also a quartic potential without double zero) . The support of the zeroes of the polynomials corresponds 
to the thick lines. There are three free real parameters ^1,^2,^3 (subject to £1 + £3 = 1, which is the 
normalization condition) . 

The third example (Fig. 20) is the case of two double 
points and double support (something that cannot happen 
for ordinary orthogonal polynomials) 

The fourth example (Fig. 19) shows that it can happen 
that not all Stokes sectors could be joined by curves T satis- 
fying the requirements for the steepest-descent method: here 
£\ € M+ is arbitrary, and pi, P2 £ are also arbitrary, but 
with the condition 5(p2) > ^s(pi)- The Stokes sector on the 
right (represented by the edge of the hexagon bordering the 
white area) cannot be joined to the other sectors because 
there is a river of h < "underwater" in between (the criti- 
cal value at the point c of intersection of the dashed lines has 
h(c) < 0). 

All these examples are admissible for our asymptotic 

Figure 18: An example of reconstruction 
study in the sense of Def. 4.2; even without a formaliza- of admissible triple of unusual topo i ogy 

tion of the rules that make an admissible triple, the reader (numerical output). 
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Figure 19: An example of admissible triple where the Stokes sectors are not all connectible by admissible 
paths (the part within the graduated square is a numerical output). Here 5(p2) > ^(pi)- 

should have no difficulty in imagining and drawing on paper 

even very complicated decorated clock-diagrams that correspond to a situation as in Def. 4.2. 

6 Conclusion 

The first two parts (Sec. 3 and 4) have shown that any admissible triple (Def. 4.2) is associated to the 
asymptotics with respect to varying weights of certain pseudo-orthogonal polynomials. Some concluding 
remarks that follow from the construction are in order; 

• there might be different orthogonal polynomials that have the same asymptotic. This can happen 
if the potential is the same (necessarily) and the asymptotic Boutroux admissible curve has two 
Stokes sectors that belong to the same connected component of {h(x) > 0}. Indeed in this case 
the jump on the contour joining them and remaining in the positive-/i region will be exponentially 
close (and uniformly) to the identity in the large N limit, thus becoming irrelevant. 

• If two Stokes sectors cannot be joined by a curve in y + U B (which happens for example in Fig 
19) then the corresponding finite-n RHP for the orthogonal polynomials cannot have a jump on a 
contour joining them. 

• In the asymptotic regime the change of Stokes' parameters (the x's) is isospectral; indeed it simply 
corresponds to choice of a different Stokes-Kirchoff 's differential ij, hence of a different twisting of 
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the line-bundle of the s. It is easy (but we don't do it here for brevity) to write explicitly the 
conjugating matrix for the spectral -problem associated to the spinors (eq. 3-55) 

The last part (Sec. 5) has shown that there can be admissible triples whose branch-cut structure B can 
be basically as complicated as one may wish and it is (topologically) a forest of loop-free trivalcnt trees; 
the zeroes of the pseudo-orthogonal polynomials then accumulate on B and we have (not completely 
rigorously) specified the asymptotic density of zeroes along these arcs in Sec. 4.6. 

6.1 Outlook 

The definition of admissibility (Dcf. 4.2) is forced upon us be- 
cause the jump matrices (Stokes' matrices) in the Riemann- 
Hilbert problem for orthogonal polynomials are all upper 
triangular. In a more general setting, e.g. as in the study 
of Painleve II equations in the asymptotic regime [19], there 
can be also lower triangular Stokes matrices: in this case we 
could relax the notion of admissibility. The key fact that 
drives the steepest-descent method is that the upper-lower 
triangular form of the jump matrix should be related to the 
signs of h on the two sides of the branch-cut (i.e. if h < 
on both sides then the jump should be upper triangular, if 
h > then the jump should be lower triangular). 

This means that one could use any Boutroux curve to 
model the asymptotics of some RHP by choosing the jump- 
matrices of the appropriate form according to the signs of h. 
This seems a promising avenue of research that we intend to 
pursue in a different publication. 

Another extension of the present setting that we set out 
to explore is the inclusion of hard-edges and more general po- 
tentials with rational derivative; while the heuristic approach 
remains unchanged (and in fact the first part of this paper is already developed to the full extent of this 
generality) several changes need to be made in the steepest descent part and in the analysis of critical 
trajectories. 

It seems that the introduction of hard-edges is the simplest generalization; this requires the use of a 
different local parametrix near each hard-edge (built out of Bessel functions) and the modification of the 
construction of admissible curves (some ingredients for this generalization can be found in [28] ) . All this 
will be dealt with in a forthcoming publication. 




Figure 20: An example of reconstruction 
of admissible triple of unusual topology 
with two double points: the constraint is 
that ^s(pi) = $s(p2) + ^s(Po)- The shaded 
area is the where h < (the "underwa- 
ter"). The dashed curves are at posi- 
tive equal height h (the dash-dot one is 
higher). Not a numerical output. 
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